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ABSTRACT 


The influence of placing interleaves between fiber reinforced plies 
in multilayered composite laminates is investigated. The geometry of 
the composite is idealized as two-dimensional, isotropic, linearly elas- 
tic media made of a damaged layer bonded between two half-planes and 
separated by thin interleaves of low extensional and shear moduli. The 
damage in the layer is taken in the form of a symmetric crack perpendicu- 
lar to the interface and may extend up to the interface. The case of an 
H-shaped crack in the form of a broken layer with delamination along the 
interface is also analyzed. The interleaves are modeled as distributed 
shear and tension springs. Fourier integral transform techniques are 
used to develop solutions in terms of singular integral equations. An 
asymptotic analysis of the integral equations based on Muskhelishvili ' s 
techniques reveals logarithmically singular axial stresses in the half- 
plane at the crack-tips for the broken layer. For the H-shaped crack, 
similar singularities are found to exist in the axial stresses at the 
interface crack-tips in the layer and the half-plane. The solution of 
the equations is found numerically for the stresses and displacements by 
using the Hadamard's concept of direct differentiation of Cauchy inte- 
grals as well as Gaussian integration techniques. 

Numerical results indicate that interleaves increase the inter- 
facial damage tolerance and significantly relieve the stresses in the 
undamaged plies. Interface cracks (H-shaped cracks) have a stable 
growth with the mode I opening stress becoming compressive after a small 
growth. Additional interface crack extension is due to shear stresses 
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(mode II) only. In order to recommend an interleaf thickness-to-layer 
width ratio, the influence of relative material properties, structural 
weight and stress reduction is studied. For example, for a typical 
Graphite/Epoxy laminate interlayered with thin interleaves made of ther- 
moplastic film, an interleaf thickness of the order of fifteen percent 
of the layer thickness is recommended. 
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CHAPTER 1 


INTRODUCTION 

Laminated fiber reinforced composite materials such as Graphite/ 
Epoxy are being used extensively in aircraft structures and are replac- 
ing many metallic components. This is mainly due to their potential for 
reducing weight and the capacity to be tailored to optimize the struc- 
tural strength and stiffness. However, due to the fiber and matrix 
interaction and multi-ply configuration of such composites, the mechani- 
cal behavior is quite complex and has challenged the designer with a new 
class of problems. One particular area which has received considerable 
attention in the past decade has been their low tolerance to interfacial 
damage. This type of damage, frequently caused by impact, is a common 
and an unavoidable occurrence during manufacturing, maintenance and ser- 
vice of aircraft structures . 

One suggested method, discussed by Masters [1] and Sun [2], to 
improve interfacial damage tolerance is to place thin films of adhesive, 
called interleaves, between those plies where delamination is more 
likely to occur. This concept is illustrated schematically by Masters 
[1, Figure 2] and is reproduced here in Figure 1. An interleaf typi- 
cally has a large shear failure strain and remains a separate layer 
between plies after curing, unlike the epoxy matrix used in the prepreg 
plies. Experimental results reported by Masters [1] and Sun [2] do 
indeed indicate that such adhesive layers are effective in reducing the 
size of interply delamination, as well as in increasing the load re- 
quired to initiate delamination. One appropriate view of this problem 
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Figure 1. Interleafed Prepreg and Cured Laminate Illustrating Retention of Discrete Layer of Tough 

Resin (Reproduced from [1] by permission of Mr. John E. Masters, American Cyanamid Company, 
Stamford, CT 06904) 
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is that impact typically causes transverse matrix cracks in the 90° 
plies due to the tensile strains caused by bending. These cracks inter- 
sect perpendicular to the 0° plies and may initiate delamination at the 
ply interfaces. If the ply thickness is small compared to the laminate 
thickness and the damaged ply is away from the mid-plane on the tensile 
side, the tensile strain gradient will be small over the ply thickness. 
The approximation of a uniform tensile strain then should give a good 
measure of the influence of the crack and the interleaf on delamination. 
An analytic solution of this problem has, however, not been given ade- 
quate treatment. A detailed solution is certainly necessary in order to 
develop an understanding of the influence of the various material para- 
meters on the behavior of interleafed composites. This is the primary 
focus of the present study. This problem could be analyzed by using 
numerical techniques such as finite element methods. But such an 
approach makes it difficult to conduct parametric studies, find stress 
singularities, and investigate crack growth behavior, which are impor- 
tant to the designer in determining the strength of the composite and 
the stress (strain) required for the initiation of damage. The proce- 
dure to be followed in the present study is based on a more complete 
solution of the linear equations of elasticity along with some appro- 
priate simplifying approximations. 

A number of related problems have been solved based on various 
forms of two-dimensional or axisymmetric approximations. Sneddon and 
Srivastav [3] solved the problem of a transverse crack in a strip of 
finite width. Gupta and Erdogan [4] considered the problem of two symme- 
tric edge cracks in an infinite strip. Hilton and Sih [5] studied the 
problem of a strip bonded to two half-planes of different materials with 
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a crack perpendicular to the interface. Bogy [6] considered the same 
geometry as in [5] to show analytically the dependence of the solution 
on material parameters. Both the solutions [5] and [6] allowed only 
embedded cracks within the strip. Ashbaugh [7] and Gupta [8] recon- 
sidered this problem with the added condition that the crack could propa- 
gate up to the interface. Erdogan and Backioglu [9] solved the fracture 
problem of a composite plate consisting of perfectly bonded parallel 
load carrying laminates and buffer strips. Gecit and Erdogan [10] 
relaxed the property of a perfect bond between plies and studied the 
effect of the thickness and elastic properties of adhesive layers 
between plies in laminated structures. The problem of two dissimilar 
elastic bonded half-planes containing a perpendicular crack terminating 
at the interface was studied by Cook and Erdogan [11]. Goree and 
Venezia [12] extended this study to include an interface crack that 
grows along the interface and crosses over the interface. Lu and 
Erdogan [13] looked at a related problem for the geometry of two dissimi- 
lar infinitely long but finite width strips . 

The intent of this study is to develop a fracture and crack growth 
model, and to investigate the influence of an interleaf on damage growth 
in laminated composites. The composite is approximated by two isotropic 
half-planes separated from a finite width layer by thin interleaves 
(Figure 2) . This is an idealization of a general laminated composite 
where one concentrates on a single damaged layer, while the outer layers 
are approximated by half-planes with average elastic properties. A uni- 
form tensile strain is assumed to be applied to the composite in the 
y-direction. A typical composite is actually multilayered and ortho- 
tropic. The approximation of isotropic, linearly elastic media gives a 
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more tractable formulation and, along with a proper choice of material 
constants and thickness ratio, should assist in understanding the influ- 
ence of the interleaves on damage growth. 

The interleaf is modeled as uncoupled distributed tension and shear 
springs. This is the most sophisticated model, other than an elastic 
continuum, available in the literature for adhesive layers. Gecit and 
Erdogan [10] used this spring model to solve the problem of periodically 
arranged dissimilar layers separated by thin adhesive layers with cracks 
perpendicular to the interface. However, this approximation was used 
only for the case of embedded cracks in the layers. For the case of 
cracks up to the interface of the layers, the spring model was replaced 
by a continuum model. This made the formulation and the solution cumber- 
some. For including delamination along the interface, which was not 
considered in their study, the continuum model for the adhesive layers 
would involve further complexity in algebraic manipulations and analysis. 
Other mathematical difficulties involved in including delamination would 
be the classical singular and oscillatory stresses near an interface 
crack-tip. The stresses undergo infinite reversals of sign as the 
crack-tip is approached, and it is also implied that the crack surfaces 
overlap near the crack-tip which is physically inadmissible. A short 
literature survey of problems with such oscillatory stresses is given by 
Cominou in reference [14]. Cominou [14] reconsidered the problem of a 
traction free interface crack between two dissimilar half-planes in a 
tension field in an attempt to explain the oscillatory stresses near the 
crack-tip. She assumed that the crack was not completely open but that 
the faces were in frictionless contact near the crack-tips. She solved 
the resulting integral equations for the length of the contact zone and 
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obtained crack-tip stresses free of oscillatory singularities. Also 
recently, Gautesen and Dundurs [15] obtained an exact solution for the 
integral equations developed in [14]. Knowles and Sternberg [16] did 
not agree with the assumptions of Cominou [14] and solved the problem by 
using the non-linear theory of elastostatic plane stress. The crack was 
found to open smoothly near the crack-tips, where the stresses were 
singular but not of oscillatory type. The spring model, on the other 
hand, does not involve such controversy and also simplifies the mathe- 
matical nature of the model. This simplification is made at the expense 
of approximating the shear and normal stresses to be constant through 
the thickness of the interleaf. An interleaf of small thickness and low 
elastic moduli as compared to those of the plies justifies such an 
approximation [10]. A detailed discussion of these approximations is 
presented later in this report. The resulting stresses in the adhesive 
are finite everywhere, while stresses with logarithmic singularities are 
shown to occur in the half-planes and the layer at the interface 
crack-tips . 

Based on these approximations, a general formulation is developed 
for plane strain and generalized plane stress . The displacement and 
stress fields are expressed in terms of Fourier transforms and, by using 
Fourier inverse transform techniques, the solution is obtained in closed 
form in terms of integral equations . These equations are solved by 
reducing them to simultaneous linear equations using Hadamard's concept 
of direct differentiation on Cauchy integrals [17] and Gaussian integra- 
tion techniques [18]. 

Three cases, depending on the extent of damage, are studied in this 
work. In the first case the center layer is assumed to have a symmetric 
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traction-free embedded crack along the x-axis; in the second case the 
crack touching the interface is analyzed, and lastly damage extending in 
the form of symmetric delaminations along the interfaces (H-shaped 
crack) is examined. 

The behavior of stresses at critical locations, for example, at the 
crack-tips, is studied to understand the influence of the relative mater- 
ial properties and geometry of the interleaf and the plies. For the 
case of a transverse embedded crack, the stresses have a square root 
singularity. For the cracks touching the interface, the axial stresses 
in the half-plane at the crack-tip are logarithmically singular. For 
the H-shaped crack, the axial stresses in the layer and the half-plane 
at the interface crack-tip are logarithmically singular. The study of 
these singular stresses and the bounded interleaf stresses are helpful 
in analyzing the influence of the relative elastic properties and geome- 
try of the interleafed composite on interface debonding, extent and sup- 
pression of delamination, and stress relieving in the undamaged plies. 

The results provide useful information that can assist the designer in 
the selection of the material and geometrical parameters of the inter- 
leaf for a particular baseline laminate. 

Certain commercial materials and their manufacturers have been men- 
tioned in this report for a practical discussion of results . The use of 
these materials in this report is not an official endorsement of these 
materials or manufacturers, either expressed or implied, by Clemson 
University. 


CHAPTER II 


FORMULATION AND SOLUTION TECHNIQUES 


Consider a laminated composite (Figure 2) in plane strain or gen- 
eralized plane stress, consisting of a single layer (Material 1) of 
width '2h', Young's modulus Ei and Poisson's ratio Vi; and half-planes 
(Material 2) having Young’s modulus E 2 , and Poisson's ratio v 2 ; sepa- 
rated by thin interleaves of thickness 't'. Young's modulus E 3 and 
Poisson's ratio v 2 . This configuration represents a multilayered 
composite with interleaves where one concentrates on a single layer with 
damage and approximates the outer undamaged layers by half-planes with 
averaged elastic properties. The damage considered is a centrally 
located crack of length '2a' along the x-axis in the layer. This crack 
may extend up to the interface (a = h) to represent a broken layer. A 
de lamination of length '2c' may grow symmetrically along the interface. 
The geometry is symmetric about the x- and y-axes. 

The composite is assumed to be loaded parallel to the y-axis with 
uniform remote stresses p t and p 2 as shown in Figure 2. The remote 
stresses are related such as to give uniform remote strains in the 
y-direction. Hence 


Pi Ei 



for generalized plane stress, and 
Pi (l-vf)Ei 

— = (l.b) 

P 2 (l-v?)E 2 


for plane strain. 
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The solution to the problem of a stress-free crack and loads applied far 
away from the crack can be obtained by superposing the solutions 
(Figure 3) of two problems. Let Oj be the solution to the problem of 
the composite without cracks and loaded with uniform remote stresses pi 
and p 2 related by equations (1). The solution of this problem is simply 
a uniform tensile strain throughout the composite. Assume Ojj to be the 
solution to the problem of the composite with damage and having no 
remote applied loads but with a uniform compression of magnitude pi on 
the transverse crack surface. 

The complete solution of the problem is hence given by 

°total = °I + °II* 

The intent of this study is to find the solution Ojj. The solution is 
found by the use of Fourier transforms and singular integral equation 
solution techniques [17, 18]. 


Stress and Displacement Field Equations 
The stress and displacement field equation for a cracked layer and 
a half-plane are developed in this section. 

The basic equations for two-dimensional linearly elastic, isotropic 
medium in the absence of body forces are given in [19]. The governing 
equilibrium equations are 


do do 

_xx _xy = o 

dx dy 


(3. a) 



dO 

- ... YY _ 
dy 


0 , 


(3.b) 


where o^, Oyy and o X y are the components of the stress tensor. The 
strain-displacement relations are 


_ 

E xx “ dx’ 


(4. a) 
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_ lY 
E yy 8y* 


3u 3v 
= * 3x* 


'xy 3y 


(4.b) 


(4.c) 


where e**, Eyy and E X y are components of the strain tensor; and u and v 
are the displacement components in the x- and y-direction, respectively. 
The stress-displacement relations are 


7“ o .... * 


K+l 


3u 3 -k 


3v 


2y xx 2 (k-1) 3x 2(k- 1) 3y’ 


7 “ 0 . = 


3-k 




K+l 3v 


2y yy 2 (k-1) 3x 2(k- 1) 3y* 


(5. a) 


(5.b) 


1_ n _ 1 ,3u 3v 

2y xy “ 2 l 3y 3x ; ’ 


(5.c) 


where y is the shear modulus, k = 3-4v for plane strain, k = (3-v)/(l+v) 
for generalized plane stress and v is Poisson's ratio. 

Using the stress-displacement expressions of equations (5), the 
equilibrium equations (3) can be expressed in terms of displacements as 


~ 2 , _2 3 / 3u , 3v. . 

U k- 1 3x ^3x 3y^ ” * 


(6. a) 


_2_ 1_ ( lu + lY 2 

k— 1 3y ^3x 3y ; 1 


(6.b) 


where 


V* = 


3 2 , 3 2 

3x 2 + 3y 2 


Elastic Layer with a Crack 

The field equations of a layer with a crack can be obtained by 
adding the expressions obtained for an uncracked layer (denoted by 
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superscript 'a') with those of an infinitely large solid with a crack 
(denoted by superscript 'b'). 

For an elastic layer without a crack, assume that the displace- 
ments are given by the following Fourier transforms, 

«• 

a 2 f 

ui(x,y) = - J F a (n,x) Cos ny dn, (7. a) 

0 


v?(x,y) = J J G a (n,x) 
0 


Sin ny dti. 


(7.b) 


where F a and G a are functions determined by the boundary conditions 
imposed in a particular problem. Substituting equations (7) in equa- 
tions (6), the mixed differential equations obtained are 


k +1 9 2 F 
a 

k- 1 3x 2 



+ 


2 

K-l 



(8. a) 


3 2 G 
a 

‘ 3x 2 




+ 




(8.b) 


Noting that the strip is of finite width and symmetric about the x- and 
y-axes, the solution of equations (8) is 

F a (n,x) = -[fi(ri) - gi(tl)] M - gi(n)x Cosh nx, (9.b) 

G a (n,x) = [f i (n) + — 1 giOO] ^skjx + gi(n)x Sinh T)x> (9. a) 

where f t and gi are two unknown functions determined by the boundary 
conditions of the problem. 

For an infinitely large elastic solid with a crack, assume the dis- 


placements to be given by the following Fourier transforms, 
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u^(x,y) = ~ | F b (S,y) Sin £x d£, 
0 


(10. a) 


v^(x,y) = ^ J G b ($,y) Cos £x d£, 
0 


(10. b) 


where Ft, and G b are functions determined by the boundary conditions 
imposed in a particular problem. Substituting equations (10) in equa- 
tions (6), the mixed differential equations obtained are 


3 2 F. 


k+1 


3G, 


r2p f — 2 _ 

3y 2 K-l D k- 1 ’ 3y ’ 


.. 3 2g k 
k+ 1 b 


3F l 


Si 17* - + Si ^ - °* 


'3y 


(11. a) 


(11. b) 


Noting that the elastic solid is symmetric about the x- and y-axes and 
that the stresses and displacements should vanish as y -*• «•, the solution 
to equations (11) is 


F b a,y) 


*i(S) 

~T~ 


, k- 1 , . -£y 

(-5 — Sy)e 


(12. a) 


G b (C,y) = 


*1(0 




$y)e 


-Ky 


(12. b) 


where $i(5) is an unknown function determined by the boundary conditions. 
The displacement field of a cracked layer is now given by 

ui(x,y) = u?(x,y) + u^(x,y), (13. a) 

Vt(x,y) = v?(x,y) + Vi(x,y). (13. b) 

Substituting equations (7) and (10) in equations (13) and using equa- 
tions (9) and (12), the displacement field of a cracked layer is written 


as 
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Kl-1 


Ui(x,y) = - J J{“ [f 1 Cti) — giCn)] Sinh(nx) 


+ xgi(n) Cosh (nx) } Cos (ny)dn 

2 f ^iCO K l“l _£y 

- - J — £ — (~2 — ■ 5y)e Sin (£x)d5, 


(14. a) 


Kj+1 


vi(x,y) = J [fi(ri) + — — gi(n)] Cosh(nx) 


+ xgi(n) Sinh (tix)} Sin (ny)dn 
2 " *i«) Ki+1 _r y 

+ - J — j? — (" 2 — + 5y)® ** Cos ($x)d£. 

0 


(14. b) 


Using the strain-displacement relations (4) and the 
stress-displacement relations (5), the corresponding stress field is 
given by 

0 1 

^; y) = | [fi(n) Cosh (nx) + nxgi(n) Sinh (nx)] Cos (ny)dn 

0 

«o 

“ J J ^i(5)(Ky)e' ?y Cos (?x)d?, (15. a) 

0 

a 1 00 

= f J { [f i(n) + 2 gl (n)] Cosh (nx) 
o 

+ nxgi(n) Sinh (nx)} Cos (ny)dn 

o» 

-Jj ♦i(C)(l%)e' 5y Cos ($x)d$, (15. b) 

0 
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XV 2y| V) = “ J {[fi(il) + gi C ti) ] Sinh (nx) 

0 

+ T)xgi(n) Cosh (tix) } Sin (ny)dn 


?y0iC?)e"^ y Sin ($x)d£. (15. c) 

0 


Elastic Half-Plane 

Assuming the same form of displacement given by equations (7), the 
solution for the elastic half-plane is given by equations (8). Noting 
the symmetry of the half-plane about the y-axis and that the stresses 
and displacement should vanish as x -+ the displacement field equa- 
tions of a half-plane obtained are 

°° 1 

u 2 (x,y) = | J {^[f 2 (h) + “ 2 — gj(n)] + Xg 2 (n)}e T1X Cos (ny)dn, (16. a) 
0 

00 £ 2 -f 1 

v 2 (x,y) = ~ J E f 2 (n ) - “ 2 — S 2 ( h ) ] + xg 2 (n)}e 1)31 Sin (tiy)dn. (16. b) 
0 

By substituting the displacement field equations (16) in the 
stress-displacement relations (5) , the corresponding stress-field equa- 
tions are given by 

a 2 °° 

= -J J [f 2 (h) + uxg 2 (n)] e _T)X Cos (ny)dn, (17. a) 

0 

c 2 00 

= J J [f 2 (n) + (t)x— 2)g 2 ( T| ) ] e _T|X Cos (Tiy)dn, 

0 


(17. b) 
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X ^ ,Y) - -J J [fi(n) + (nx-l)g 2 (n)] e“ nx Sin (ny)dn. (17. c) 

0 

The unknown functions f tl g 1( $ lt f 2 and g 2 in the displacement and 
stress field equations (14-17) are determined by applying appropriate 
boundary and continuity conditions relating to the problem. The solu- 
tion is obtained in terms of singular integral equations which can be 
solved by standard numerical techniques [17, 18]. 

Derivation and Solution of the Integral Equations 

Three cases related to the extent of the damage are investigated. 

The first case deals with a symmetric transverse crack of length '2a' 

(a < h) located centrally along the x-axis. In the second case, the 
crack can extend upto the interface (a * h) and represents a broken 
layer. The last case accounts for a crack up to the interface (a = h) 
along with splitting of length '2c' located symmetrically parallel to 
the y-axis at the interface (x = ±h) (Figure 2) and is more commonly 
called the H-shaped crack. The nature of the damage effects the continu- 
ity and boundary conditions of the problem. Hence, each case is studied 
separately. The solution is derived in terms of singular integral equa- 
tions. An asymptotic analysis is carried out to study the singular 
behavior of stresses and displacement slopes in the vicinity of the 
crack-tips. Numerical techniques are developed to solve the integral 
equations for the unknown functions. 

Crack within the Layer 

As discussed in Chapter 1, the interleaves are modeled as 
uncoupled distributed tension and shear springs and the stiffness of 


these springs is given in [10] as 


Eo 


t"» 

(18. a) 

v 3 



(18. b) 


where 

E 3 

Eg a , 2 

0 l-v 3 

for generalized plane stress, 

EaCl-Va) 

“ (l+Va)(l-2va) 
for plane strain. 

By modeling the interleaf as springs, the transverse and shear 
stresses are constant through the thickness of the interleaf. The con- 


tinuity conditions along x 

= h are then written as 


* °xy (h ' y) ' 

o < 1 y 1 


(19. a) 

°xi (h > y) =0 « (h - y) - 

0 < | y | 

< 

(19. b) 

°xx (h,y) = K n f U2(h,y) 

- Ui (h,y) ] , 

o < 1 y 1 < «, 

(20. a) 

° X y(h,y) = K s [va(h,y) 

- Vi (h,y) ] , 

o < 1 y 1 < 

(20. b) 

The homogeneous boundary conditions along 

y = 0 are 


°xy( x, 0) = 0> 

| x | < h, 


(21. a) 

°xy (x,0) = °» 

| x | > h, 


(21. b) 

v 2 (x,0) = 0, 

| x | > h. 


(21. c) 


The mixed boundary conditions along y = 0 are 
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Vi(x,0) =0, a < | x | < h. (23) 

The crack surface traction is a constant for this study and is given by 
p(x) = pj as defined by the solution Ojj in equation (2). However, the 
problem is formulated here for the general case of crack surface trac- 
tion as a symmetric function of 'x' , where p(x) * pi is a special case. 

Substituting the stress and displacement field equations (14-17) 
in the continuity conditions (19) and (20), and taking the inverse 
Fourier transforms in 'y' yields the following simultaneous equations. 


aii fi(n) + ai 2 (n) gi(ti) + a^ f 2 (ri) + a^ g 2 (n) = Di(ri), 

(i-l,...,4), (24) 


where aij, (i,j = 1,...,4) and Di, (i = 1,...,4) are given in Appendix B 
and the integral formulas in Appendix A are used. Equation (24) can be 
solved simultaneously for fi and g x in terms of the unknown functions 
Di, which are infinite integrals of the function $i(£) to obtain 


fi(n) 


FiOj 2 
® 1 1® 2 2 


F 2®1 2 
<*21012 



(25. a) 


gi(n) 


FjOi i 
011022 


Fl ° 21 -nh 
e , 

«2 l«l 2 


(25. b) 


where ai j , (i = 1,2, j = 1,2), and Fi, (i = 1,2) in terms of Di, 

(i = 1,...,4), are given in Appendix B. 

The homogeneous boundary conditions (21) are identically satisfied 
by the stress and displacement field equations (14-17). 

The mixed boundary condition (23) may be expressed as 


v(x) = Vi(x, 0 ) = 



0 


♦ i($) Kt+1 . . 

—z (— 2 — ) Cos = a<|x|<h. (26) 
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Taking the Fourier cosine transform of equation (26) gives 

♦i($) 2 f 

J v(t) Cos$t dt. (27) 

0 

Substituting equation (27) into the expressions for Dj, (i = 1,...,4), 
given in Appendix B, gives 
a 

( Kl +1)D. = j H^Ti.t) v(t ) e “ ,l(h_t) dt, (i = 1, . . • ,4) , (28) 

-a 

where H^C^t), (i = are given in Appendix B. Using these 

expressions for D^, the functions F^, (i = 1,2), used in equation (25) 
can be expressed as 
a 

F i ( n) = j E.(n,t)v(t)e" T,(h_t) dt, (i = 1,2). (29) 

-a 

The functions E^, (i = 1,2), are given in Appendix B. 

The mixed boundary condition (22) may now be written as 

CCx.O) - ~ 

= lim - - [ *i(£)(l+£y) e Cos £x d$ 

2Ui , u J 

y-*-0 0 


+ 


2 


•© 

J 


{(fi(ri) + 2gi(Ti)] Cosh nx + nxg!(n) Sinh iix} dn 


0 


p(x) 

~ " 2yj 5 



(30) 


With defined by equation (27), the first integral in equation (30) 


becomes 
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lim - “ f $i(£)(l+Sy) e Cos £x d$ 
y*o* J Q 


a 


= ttCkx+I) J vCt) [ (t-x) 2 + (t+x) z] dt ‘ 


(31) 


The even symmetry property of the crack opening displacement, v(x) 
simplifies the right hand side of equation (31) to 
a 

tKkx+I) J v(t) f (t-x) 2 + (t+x) a] dt 


= 2 f .v(t) . 

TrCKi+l) J (t-x ) 2 


(32) 


-a 


Substituting equation (29) in equation (25), the unknown functions 
and g i can be algebraically written as a function of the displacement 
function v(x) • Substituting these results in the second integral of 
equation (30), and equation (32) in the first integral, a singular 
integral equation in terms of the unknown displacement function v(t) is 
obtained as 


I (t-x)* dt + J v ( t ) dt - ~ 


ir(l+K 1 )p(x) 

4Vi 


, | x | < a, (33) 


-a 


-a 


where 


Kn(x,t) = | k 11 (t,x,tj) e 

0 


The expression for k 11 (t,x,n) is 


-n(h-t) 


dn. 


(34) 
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kn(t,x,n) = {[o22”2a 2l ) Cosh nx - tix o 2l Sinh T)x] E x 

+ [(-Oi2+2on) Cosh nx + nx a 2 j Sinh nx] E 2 }e 


-nh 


/(otl l« 22 -«lll* 12 ). (35) 

where <*ij» (i = 1,2, j = 1,2), and E^, (i = 1,2), are given in Appendix 
B. 

Using integration by parts, the singular integral equation (33) can 
be expressed in terms of the slope of the crack opening displacement, 
defined as 


G(x) = 


_ 3v(x) 


9x 


(36) 


The integral equation then has the following form 


v(a) F(x) + J dt + J G(t) K ll (t,x)dt = - 

-a -a 


irp(x)(l+K 1 ) 

i 


x I < a. 


where 


(37) 


F(x) = J [{fn [(ot 2 2 “ 2 o 2 i) Cosh nx • nx Bn Sinh nx] 

0 

+ f 2 al(”®i 2 + 2on) Cosh nx + i)x a M Sinh nx]}e 


-nh 


/ («ii <*22 - O 21 O 12 ) “2 Cosh nx e dn, 


(38) 


and fit and f 2 * are given in Appendix B. 
The integrand £n(t,x) is 


Kn(t,x) * J kn(t,x,n) dn, 


(39) 
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where 


9E! 

kn(t,x,n) = {[(a 22 -2a 21 )Cosh nx - T)x o 21 Sinh T)x] 


3E ; 


2 -nh 

+ [(‘Ou + 2on) Cosh nx + HX o u Sinh tjx] j£“}e 


/ ( a ll a 22 ” 821812) , ( 40 ) 

and aij, (i = 1,2, j = 1,2), and E^, (i = 1,2) are given in Appendix B. 

Since the vertical displacements, v(±a) at the crack tips must be 
identically equal to zero, the integral equation (37) reduces to 

f G(t) f - -irp(x)(l+ Kl ) 

J + J G(t) Kn(t,x)dt = — , | x | < a. (41) 


-a 


-a 


All the special cases can be recovered from equation (41). For 
example, for the interleaf thickness t = 0, equation (41) is identical 
to that given by equation (17) in [8]. The problem of a layer, with a 
crack situated perpendicular to the interface, bonded between two 
half-planes of elastic properties different from the layer was solved by 
Gupta in [8]. Equation (17) in [8] can also be obtained by making the 
interleaf infinitely stiff, that is, make the elastic moduli y 3 and E a 
approach infinity. Note that the interleaf can be made stiffer by 
either increasing its elastic moduli or by decreasing the interleaf 
thickness. For the other special case of letting the interleaf thick- 
ness approach infinity, equation (41) is identical to equation (19) in 
[4] . The problem of an infinite strip with an internal crack perpendic- 
ular to the stress-free longitudinal edges was solved as a special case 
by Gupta and Erdogan in [4] . 

It can be shown that kxi(t,x,Ti) is a bounded and exponentially 
decaying function. Hence Kn(t,x) is also bounded and the singular 
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integral equation is of the Cauchy-type. The singular behavior of the 
slope function, G(t) may be analyzed by studying the dominant part of 
the singular integral equation (41), which can be expressed as 
a 

1 SP dt = I x | < a, (42) 

-a 


where B(x) contains all the bounded terms of equation (41). 

Following the technique explained in [8], G(t) is assumed to have 
an integrable power singularity at t = ±a which can be written as 


G(t) = 


H(t ) 

(a 2 -t 2 ) T 


H(t)e 

(t-a)*(t+a) y 


I t | 


a. 


(43) 


where 0 < ReY < 1 and H(t) satisfies the Holder conditions [20] in the 
closed interval [-a, a]. A function f(x) is said to be Holder continuous 
in [-a, a], if for any two points Xi, x 2 in [-a, a] | f(x 2 ) - f(xx) | < 

A | x 2 -Xi f is satisfied. A and p are positive constants and 0 < p < 1, 
'A' is called the Holder constant and 'p' is termed as Holder index. 

Consider the sectional ly holomorphic function 


♦ <z) -if pi dt - i f 

ii J t-z ir J 


H(t)e 


niY 


-a 


-a 


(t-a)*(t+a)\t-z) 


dt. 


(44) 


According to ([20], Chapter 4) 

»iY 


. H(-a) e 

«(z) = — ± r 


, 0 .Y SirnrY , . ,Y 
(2a) (z+a) 


1 _ Hfal _X 


.Y SintY , . 

(2a) (z-a) 


y + 0 o (z) . (45) 


where $ 0 (z) is bounded everywhere except at the end points z = ±a, where 

(46) 


H 0 (±a) 

Uo(z)l<; Re(Yo) < Re(Y). 


I z±a f 9 


Following [8] and using the Plemelj formulas, 
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/(*) - T(x) = 2iG(x) , | x | < a, (47. a) 

a 

♦ + (X)+#"(X)-J J ^ dt, I x | < a, (47. b) 

-a 

, f s _ H(-a~)Cot ttY H(a1 Cot 

r(.XJ y r y 

(2a) (a+x) fl (2a) a (a-: 

* 

where $o(x) is a bounded function and has behavior similar to (z) . 
Substituting equation (48) in equation (42) and letting x -*■ a or -a, 

Cot ttY = 0 (49) 

is obtained. Equation (49) is called the characteristic equation and 
in the acceptable range of 0 < Re Y < 1 has the root 

Y * (50) 

Using the above value of Y, equation (43) gives the slope of the 
displacement as 

G(t) = | t | < a. (51) 

✓a 2 -t 2 

Substituting equation (51) in equation (41), the singular integral equa- 
tion may then be expressed as 

X I < a. (52) 

-a 

Since the index of the integral equation (41) is +1 [18], the solution 
will contain an arbitrary constant, which is determined by the single 
valuedness condition 
a 

J G(t)dt = 0. (53) 

-a 


f It zz + Kn(t,x)]dt = 

J vV-t 2 1 x 


*(1+Ki)p(x) 
4 \ii 


itY 


+ *o(x), 


a. 


(48) 
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Normalizing the variables with respect to the half crack length 


v. 


s = 


(54. a) 


r = 


a’ 


(54. b) 


G(t) = G(as) = $(s), 
p(x) = p(as) = S(r), 

H(t) = H(as) = T(s), 
equation (41) may be written as 
1 

/ + a Kn(as,ar)]ds = 

-1 


tS(r)(l+K 1 ) 

4px 


r < 1, 


(54. c) 
(54. d) 
(54. e) 


(55) 


Using equation (51), the non-dimensional slope function, 0(s) can be 
expressed as 


*(s) = , | s | < 1. 


a/l-s 2 


(56) 


The integral equation (55) is approximated by using Gauss-Chebyshev 
integration formulae and techniques [18] to give 


N 


N X, +aK " (a V al j’l 

1=1 i y J 


irS(r j )(l+K 1 ) 

4yi 


, (j = 1,2, .. . ,N-1) , 


(57) 


where 

s i= Cos , (i = 1,2,...,N), (58) 

r Cos [J 1 ], (j = 1,2, .. . ,N-1) . 


(59) 
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Using equation (56), the single valuedness condition (53) can be 
approximated as 
N 

l Y(s.) = 0. (60) 

i=l 

Equations (57) and (60) represent a N x N system of simultaneous linear 
equations which gives numerical values for i'(s) at the discrete 
points, s^ given by equation (58). 

Once ns) is determined numerically, the unknown functions fj, 
gx, f 2 » g 2 , and 0i can be evaluated using equation (24). This implies 
that the stresses and displacements can then be obtained at any point. 
The normal cleavage stress, (a + ,0) at the crack-tip, has a 

square root singularity (see Appendix C). A stress intensity factor, Kx 
is defined as 


Ki = lim /2(x-a) o‘ (x,0). 
x->a + 


(61) 


Equation (55) can be re-written (Appendix C) in non-dimensional terms as 


Kx 


U+ki) / a 


(62) 


In the numerical solution, the value of the non-dimensional func- 
tion T(s) can be determined only at discrete points given by equation 
(58) and the end points (±1) cannot be chosen as one of these discrete 
points. Krenk [21] developed a summation formula by using properties of 
Chebyshev polynomials to give the extrapolated value of T(l) as 


N 


f(l) - n I 

i=l 


Sin[^ (2i— 1 )tt ] 

Sinflg 1 it] 


t(s.). 


(63) 
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The routines used to solve the kernel Kix(x,t) and the discussion 
on the convergence of the results is given in Appendix E. Convergence 
of the results for similar elastostatic problems obtained by using the 
techniques discussed for the present case is discussed in detail in 
[ 22 ]. 

Crack Up to the Interface 

The boundary and continuity conditions for the crack extending up 
to the interface remain the same as for the case of the crack within the 
layer as discussed in the previous section. Hence, the governing singu- 
lar integral equation (37) is valid for this case too. However, 
the Fredholm Kernel K lx (t,x) becomes unbounded for certain values of x 
and t. An asymptotic analysis reveals that K u (t,x) is unbounded as 
x -*■ h and t -*■ h or x -h and t -*■ h, simultaneously. Since the 
integrand ku(t,x,ri) in equation (37) is bounded and continuous in the 
semi-infinite interval 0 < n < », the unbounded terms in K xl (x,t) are 
due to the asymptotic behavior of kn(t,x,i]) as ii ■+ ». The part 
contributing the unbounded term may be expressed as 


K7i(t,x) = J k7i (t,x,»i)e” T, ^ 2h " t ^ dn, 
0 


with 

k7i(t,x,n) = [-4h(h-t)u 2 + n{6(h-t) + Bxh(h-t) + 2h} 

+ {[B 2 h-4+B 3 h(h-t) + B* (h-t) } ] Cosh nx 
+ [4x(h-t)T) J + n{-BxX(h-t)-2x} 

+ f-B 2 x-B 3 x(h-t)}] Sinh tix, 


(64) 


( 65 ) 


where the functions B^, (i = 1,...,4), are given in Appendix B. 
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Using the following integrals [23] 


I 


m -TjfZh-t')-. . d , x 

n e 'Sink m dn = — l (2h _ t) »_ x .l. 

at 


i 


m -n(2h-t)_ , . d“ r 

T 1 e Cosh nx dn = [ 


dt 


2h-t 

2 2 J > 


m 1 (2h-t) z -x 


the singular kernel Ku(t,x) may be written as 


KTi(t.x) = \ {-4(h-x) 2 + (12(h-x) + Bi (h-x) 2 ] £ 


dx 


+ [-2- (B i+B 2 -Bi, ) (h-x) + B 3 (h-x) 2 ] } t _ ( 2h_x ) 


(66. a) 


(66. b) 


+ \ {-4(h+x) 2 + f-12(h+x) - Bi (h+x) 2 ] } ^ 

+ [-2-(B l +B 2 -B,)(h+x) + B 3 (h+x) 2 ]} t _ ( 2 h +x) ‘ (67) 

The coefficient of v(h) in equation (37) also becomes unbounded as 
x -*• ±h. The asymptotic behavior of F(x) shows that the part contri- 
buting the unboundedness may be expressed as 


F~(x) = J f“(n,x) e' 1111 dt| , 
0 


( 68 ) 


where 


f (T),x) = n+ij Cosh ' nx ’ 


(69) 


and B^, (i = 5,6), are given in Appendix B. 

The behavior of the slope function, G(t) at the end points can now 

be found by considering the dominant part of the singular integral 
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equation which can be expressed as 

~ h 

v(h) | f“(n,x)e _T,h di| + J G(t) + KTi(t,x)]dt = B(x), | x | < h, (70) 

0 -h 


where B(x) is a bounded function. Assume G(t) to have an integrable 
singularity at t = ±h given by 


G(t) = 


HCt) 

(h 2 -t 2 )* 


H(t)e iriy 

(t-h)*(t+h)* 


t | < h, 


(71) 


where 0 < Re Y < 1 and H(t) satisfies the Holder condition in the closed 
interval [-h,h]. Consider the sectionally holomorphic function 


* (2) * i / S? dt - 

-h 


(72) 


Following the procedure outlined in the previous section and applying 
the Plemelj formulas (47) to equation (72), 


#(X) = + »t(x), I X I < h, 

(2h) (h+x) 4 (2h)*(h-x) 4 


(73. a) 


$(2h+x) = - 


H(,h ) 


(2h)*Sin irY (h+x) 


+ *?(x), h < 2h + x < 3h, (73. b) 


*(2h-x) = - 


HQ) 


(2h) y Sin ii Y (h-x) J 


1 + *?(x), h < 2h - x < 3h, (73. c) 


where $i(x), (i = 1,2,3), has behavior similar to $ 0 (z) given by 
equation (46). 

Substituting equations (65) and (68) in equation (70) and using 
equations (69) and (73) to separate the unbounded terms, the following 
equation is obtained 
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v(h) [B s e Bs ^ h X ^log(h-x) + B s e B * ^ h+X ^log(h+x) ] 

+ y 1 ' 1 ^ r (" Cos + i { 4jr (* +1 ) ” 12y + 2}] 

(2h) Sin uJ(h+x) 

+ x ^ if t" Cos v * + 9 {^TC^+1) - 12* + 2}] = B(x), (74) 

(2h) Sin *2T(h-x) 4 Z 

where B(x) is a bounded function. 

X 

Multiplying equation (74) by (h+x) and substituting x = -h and 

y 

then multiplying by (h-x) and substituting x = +h, the resulting 
characteristic equation is 

-Cos ir* + 2Z* - 43T + 1 = 0. (75) 

This characteristic equation is independent of material constants and is 
identical to the characteristic equation obtained for the solution of a 
free end problem of a strip under self-equilibrating residual stresses 

19 ]. 

There are no roots of equation (75) in the acceptable range of 
0 < Re Jf < 1, which implies that G(t) does not have a power singularity. 
Hence G(t) either has a logarithmic singularity or is bounded at the end 
points t = ±h. 

If G(t) is a bounded function, then 

G(t) = H(t), t < | h |, (76) 

where H(t) satisfies the Holder conditions in the closed internal [~h,h] . 
Again consider the following sectionally holomorphic function 

- 1 ■ i s? dt - 

-h -h 


( 77 ) 
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According to [20, Chapter 4] 

*(z) = H(h) log(z-h)-HC-h) log(z+h) + * 0 (z), (78) 

where $ 0 (z) is a bounded function tending to a definite limit when 
z -*■ ±h along any path. Following [9], using the Plemelj formulas (47), 
and noting that G(t) is an odd function, the following results are 
obtained, 

0(x) = H(h) log(h 2 -x 2 ) + *?(x), | x | < h, (79. a) 

*(2h+x) = H(h) log(h+x) + , h < 2h + x < 3h, (79. b) 

$(2h-x) = H(h) log(h-x) + tfttx), h < 2h - x < 3h. (79. c) 

ic 

The functions $i(x), (i = 1,2,3), are bounded functions. 

Substituting equation (79) in equation (70) and separating the 
unbounded terms 

v(h)[B s e Bs(h " x) log(h-x) + B s e Bs(h+x) log(h+x)] =B(x), (80) 

where B(x) is a bounded function. The above equality is possible only if 
the displacement at the end points, v(±h) is zero. This is not the case 
because the interleaf is modeled as distributed tension and shear springs, 
which allows the broken layer to displace at the end points under a uni- 
form pressure, pi on the cracked surface. Hence, G(t) cannot be bounded 
at the end points. This leaves only the possibility of a logarithmic 
singularity for G(t) at the end points. 

Assume G(t) has an integrable logarithmic singularity at t = ±h 
expressed as 

G(t) = H(t) log(j^), |t|<h, (81) 


where H(t) satisfies the Holder condition in the closed interval [-h,h] . 
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Note that since G(t) is an odd function, H(t) must be an even 
that is, H(t) = H(-t). 

Consider the sectionally holomorphic functions 
h 

0(z) = J 

-h -h 


s^t- 


t“Z 


-J 


t-z 


dt 


= 0l(z) - 0 2 (z), 


where 


0l(z) 


I 


H(t)log(h+t) 

t-z ’ 


0 2 (z) 


h 


I 


H(t)log(h-t) 

t-z ’ 


According to [20, Chapter 4], near z = h. 


Mz) = H(h) log(2h) log(z-h) + 0 o (z), 


where 0 o (z) is a bounded function tending to a definite limit 
Near z = -h, 

h 


-h 


Consider 


dt + bounded function. 


h 


fl(z) 


1 

2ni 


! 


]og(h+ t ) dt 
t-z 


-h 

Using the Plemelj formulas (47), 


function, 


(82) 
(83. a) 
(83. b) 

(84) 

at z = h. 

(85) 

( 86 ) 


ft + (t 0 ) ~ (t 0 ) = log(t 0 +h) for t 0 t [-h,h]. 


(87) 
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Also 


u + (t 0 ) “ »“(t 0 ) = -log(t 0 +h) for t 0 £ [-h,h], (88) 

if 

(log(z+h)} 2 

w(z) ■ 2 f 2 itI log(z+h)]. (89) 

Adding equations (87) and (88), 

(8+w) + - (8+«)" = 0, (90) 

and therefore the function [8(z) + w(z)] is holomorphic in the neighbor- 
hood of z * -h, that is, near z = -h, 

, {log(z+h)} * 

Q(z) * **2 [ log(z+h)] + holomorphic function. (91) 

From equations (82), (83), (84) and (91), the behavior of <h(z) near the 
end points z * ±h is written as 

♦ i (z) = H(h) log(2h) log(z-h) - [{log(z+h)} 2 - 2iri log(z+h)] 

+ holomorphic function. (92) 

Similarly the behavior of $ 2 (z) near the end points z = ±h can be 
expressed as 

* 2 (z) = H(h) log(2h) log(z+h) + [log(z-h)] 2 < 93 ) 

+ holomorphic function. 

Applying the Plemelj formulas (47) in equations (92) and (93), equa 
tion (72) is reduced, for z = x, 2h + x, and 2h - x, to give 

<Kx) = H(h) log (2h) log(h 2 -x 2 ) - [{log(h-x)} 2 


+ {log(h+x)) 2 ] + *?(x), | x | < h. 


(94. a) 
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0(2h+x) = [log(h+x)]* + H(h) log(2h) log(h+x) 

+ #*(x) , h < 2h + x < 3h, (94. b) 

0(2h-x) = (log(h-x) ] 2 + H(h) log(2h) log(h-x) 

+ $*(x) , h < 2h - x < 3h, (94. c) 

where $*(x), (i * 1,2,3), are bounded functions. 

Substituting equations (94) in equation (70), and separating the 
unbounded terms 

[-B s e B *^ h_X ^ log(h-x) - B b e B ‘^ h+X ^ log(h+x) ] v(h) 

- 4H(h) [log(h-x) + log(h+x)] = B(x), (95) 

where B(x) contains all the bounded terms. For the boundedness of the 
left hand side of equation (95), it is required that 
B s K!+l y 3 

H(h) = - — v(h) = - — — v(h) . (96) 

Noting that 

o‘ (h,0+) = o* ( h,0+) = - — v(h) , (97) 

xy xy *- 

and 

H(h) = °xy (h ’ 0+)> (98) 

then 

K i +1 h+v 

lim G(x) = — — 0* (h,0+) lim lo 8(hI^)- 

x-*±h“ Ul y x-+±h" 

This proves the assumption that the slope, G(x) of the crack opening 
displacement, v(x) function has a logarithmic singularity at the enc * 
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points is correct. Note that the coefficient of the logarithmic term is 
dependent on the shear stress at the end point, (x = h, y ■ 0+) . 

The problem of a crack in the layer extending up to the interface 
can also be viewed as two semi-infinite strips with uniform pressure, p t 
applied on the transverse axis (y = 0) and unknown finite tractions on 
the longitudinal sides (x = ±h) . These tractions depend on the elastic 
properties and geometrical parameters of the interleaf and the 
half-plane. This geometry can be compared with the known solution for a 
semi-infinite strip (obtained from the free-body diagram of the classi- 
cal half-plane problem with uniform pressure over a finite range, Figure 
4) having the following boundary conditions 


O X y( x >0) 38 

0, 



-h < x < h, 

(100. a) 

0yy(x,O) = 

-Pl» 



-h < x < h, 

(100. b) 

°xx (ih ’ y ) * 

Pi 

n 

[§ - Tan- 1 

_ 2hy . 

2h 4h 2 +y 2 J ’ 

0 < y < «, 

(100. c) 

° xy 0>.y) = 

Pi 

V 

4h* 

4h 2 +y 2 * 


0 < y < «», 

(100. d) 

V h ' y) - 

pi 

* 

4h 2 

4h 2 +y 2 * 


0 < y < -. 

(100. e) 


The solution to the above boundary value problem is given in [ 19 , 
Chapter 5] by using the Airy stress function method. The slope of the 
vertical displacement is given by 


3v(x.(n 

dx 


_ 2Pi i /h±x. 
*E log ( h-x ) ’ 


or 


(101. a) 


av(x, P) _ 2 

dx ~ E 



(h,0+) log( 


h+x 

h-x 


)• 


(101. b) 




Figure 4. Classical Half-Plane Problem with Uniform Pressure Over a Finite Range 
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It is interesting to note that both the problems have a logarithmic 
singularity at the end points x = ±h. In fact, equations (99) and 
(101. b) are identical expressions for the slope function, G(x) at x = ±h. 
A finite non-zero shear stress at the point (x = ±h, y = 0+) is the rea- 
son for the logarithmically singular slope functions in the two problems. 

The shear stress in the half-plane is discontinuous, non-zero and 
finite at x = h, y = 0. Hence, the axial normal stress is logarithmi- 
cally singular [19, 24] at x = h, y = 0 in the half-plane and is given 
by 

lim o* (x,0) = K lim log(x-h) + Order (1), (102) 

x-h + 77 x-+h + 


where K is a modified stress intensity factor given by 

a 4y 3 v(h) 

K = ~ o 1 (h,0+) . 

ir xy ut 


(103) 


The same expression is obtained by an asymptotic analysis (see 
Appendix D) . It is interesting to note that a discontinuity in the nor- 
mal loads on the half-plane does not produce unbounded stresses . 

The numerical techniques enumerated in [18] to solve singular inte- 
integral equations with Cauchy kernels are based on the exact Cauchy 
gral equations with Cauchy kernels are based on the exact Cauchy prin- 
cipal value integral expression [25] of orthogonal Jacobi polynomials, 
P n (t) given by 


} (l-t)“(l+t)M a,e) (t) 

I 


= ir Cot(air ) ( l-x) a ( l+x)° P (a,e !x) 

n 

F(n+1; -n-a-8; 1-a; “)(a>-l, 


. 2 a . +g r(«)r(n+p+i) 

T (n+a+8+1) 

B>- 1 , 01*0,1,2, . . .), | x | < 


1 , 


(104) 
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where T(a) is the Gamma function and F(a;f5;JT;z) is a hypergeometric 
series. This relationship forms the basis for solving singular integral 
equations with a simple Cauchy kernel, where the unknown slope function 
has a power singularity of the order of (a) and (&) at the end points. 
But in the present case the slope function has a logarithmic singularity 
at the end points for which a similar numerical technique could not be 
found in the literature. Hence solving equation (33) for the unknown 
displacement function, v(x) is suggested. 

Normalizing the following variables with respect to the half-layer 
width 'h', 


t 

s = h’ 

(105. a) 

X 

r = h’ 

(105. b) 

V(s) = h v(t). 

(105. c) 

Kn(r,s) = h 2 K lx (t,x), 

(105.d) 

S(r) = p(x) , 

(105. e) 


equation (33) can be re-written in the non-dimensional form as 
1 1 

J ( g.rjz ds + / h * V(s) K lx (hs,hr)ds 

-1 -1 

TfS(r)(l+ic 1 ) 

^ , I r I < 1. (106) 

Equation (106) has z ~ ry type integrands and are called strong 

1st ) 

singularities. Such integrands are classically non- integr able and 
cannot be defined even in the principal value sense. A concept used by 


40 


Hadamard [17] interprets improper integrals of such singularities in the 
finite part sense. Kaya and Erdogan [26] used this concept to solve 
elastodynamic problems such as the problem of an edge crack perpendic- 
ular to the free boundary of an elastic half-plane. 

Using the Cauchy principal value integrals [27] 


. P (t) 

J ~: x dt = -2 Q n (x)> I * I < i. 

-1 


( 107 ) 


where P n (t) and Q n (x) are Legendre polynomials of the first and the 
second kind, respectively. 

The Hadamard' s concept uses direct differentiation on Cauchy inte- 
grals to obtain finite part integrals. This is given by 



plat = f 

t-X J 

-1 


-fit} 

(t-X ) 5 


dt. 


( 108 ) 


Using the above direct differentiation in equation (107), 


P (t) 

/ (fep dt = '(lSH [x V x) ' Vl (x)l> I x I < l * (109) 

-1 


and assuming 


V( S ) * I A P (.), 

n=0 


( 110 ) 


equation (106) can be written as 

N 1 

l V " 2 i- 7 ^ {t Q n (r) - Q n+1 (r)} +j P n (s) K x x (hs ,hr)ds] 


-1 


irS(r)(l+Ki) 

4y x 


r < 1. 


(Ill) 
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A collocation method is used to solve for unknown coefficients, A n . 
The choice of collocation points, symmetrically distributed on the inter- 
val (-1,1), is not restricted but more points concentrated near the end 
points improves the rate of convergence . The roots of the Legendre poly- 
nomials are a suitable choice of collocation points to give 


N 


l A [ ^(n+l) { r , Q ( r .) - Q 
L n n l 1-r. 1 i x n i n+1 

n=0 i 


+ J 


P n (s) K 11 (hs,hr i )ds] 


-1 


irS(r i )(1+K!) 
4vi 


(i = 0,1,2, ...,N-1,N), 


(112) 


where r^ is the i** 1 zero of the Legendre polynomial Pn+i(x), and 

P N (X) = ~N~ 

21! dx 


(113) 


Equation (112) represents a [(N+1) x (N+1)] system of linear equa- 
tions which can be used to numerically evaluate the coefficients, A n . 
Using equation (110) the crack opening displacement can be determined. 
Stresses and displacements can be evaluated by finding functions f ly g x , 
f 2 , g 2 and via equation (24). 

The modified stress intensity factor, K defined by equation (103) 


is expressed numerically as 


K = - 


4y 3 h 

ut 


N 

X 

n=0 


A 


n 


(114) 


The routines used to evaluate the kernel K u and the finite inte- 
gral over [-1,1] in Equation (111), as well as the convergence of the 
results are discussed in Appendix E. 
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H- Shaped Crack 

The damage considered in this case is a broken layer (a = h) with 
delamination of length ' 2 c' along the interface (Figure 2 ). By 
including delamination, only the continuity conditions along x = h are 
different from the previous two sections, and these are given by 


= 0 ix Ch - y) > 

o < 1 y I < -, 

( 115 . a) 

o* T (h>y) “ o* o>.y). 

o < I y I < — , 

( 115 . b) 

xy xy 



°ix (h ’ y) = °’ 

0 < | y | < c. 


= K n [u*(h,y)-u 1 (h,y)] , 

c < | y | < •*. 

( 116 ) 

o* y (h,y) * 0, 

0 < 1 y 1 < c. 


* K s [v 2 (h,y)-vi(h,y)] , 

c < | y | < -. 

( 117 ) 

Defining two unknown functions, $ 2 (y) 

and $3 (y) by 


2Vi*j(y) 



k !+ i = 0 ix< h '^ - Mu>o>>y> 

- Ui (h,y) ] , 0 < | y | < — , 

( 118 ) 

2Vi^ 3 (y) 



- Kl +1 * °iy th - y) ‘ ^ |V!<h ' y) 

- vi(h,y)] , 0 < | y | < -, 

( 119 ) 

the continuity conditions ( 116 ) and ( 117 ) can be re-written as 


°ix (h * y) = °» 

0 < | y | < c, 

(120) 

2yi0 2 (y) 

1 1 


Kl+ 1 = °ix (h > y) “ K n [u ’ (h ’ y) 

- Ui (h,y) ] , | y | < c, 

(121) 

* 2 (y) = o. 

1 y 1 > c, 

(122) 

o* (h.y) - o, 

0 < | y | < c, 

(123) 
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2yiMy) . 

— = °iy (h,y) " K stMh,y) - v 1 (h,y)] > I y I < c, (124) 

My) =0, | y | > c. (125) 

Substituting the stress and displacement equations (14-17) in equations 
(115), (120), (121), (123) and (124) and taking the inverse Fourier 
transforms of the resulting equations in 'y' > and using the integral 
formulas given in Appendix A, the following equations are obtained. 

a iifi(n) + a i2gi(n) + a i 3 f 2 (Ti) + ai*g 2 (n) * L^n), 


(i = 1, • • • ,4) , 

where a^ j , (i * 1,...,4, j = 1,...,4), and L^, (i 
in Appendix B. 

Using the mixed boundary condition 
v(t) Cos(t dt. 

0 

for a s h, and the continuity conditions 


MS) 2 t 

S ~ Kl+1 J 


(126) 

1, . . . ,4) , are given 


(27) 


My) =0, I y I > C, (122) 

My) = 0, I y I > c, (125) 

the right hand side, L^, (i = 1,...,4), of equation (126) is expressed 
as finite integrals in terms of the unknown functions v(t), <t> 2 ( y) and 
My) and are given in Appendix B. 

Equation (126) can be solved simultaneously for f v and gj in terms 
of the three unknown functions v(t), $ 2 (y) and $ 3 (y) to obtain 


fl(Tl) 


Mi<x 2 2 ”M 2 a t 2 
°1 1°2 2” a 2 1®1 2 



gi(n) 


M 2 0! 1-Mi0t 2 x 


a l l a 2 2 -a 2 l&l 2 



(127) 
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where o^j, (i = 1,2, j = 1,2), are given in Appendix B, and 

h c 

Mi(rt) = [ J Ei e T *(h~ t ) v (t)dt - X 3 n | $ 2 (y) Cos ny dy 

-h 0 

c 

— X* n J $a(y) Sin ny dy] , 

0 

h c 

M 2 (n) = [ J E 2 e" 11 ^ h “ t - ) v(t)dt + X s n | $ 2 (y) Cos ny dy 

-h 0 


c 

- x s H J * 3 (y) Sin ny dy] . (128) 

0 

where E^, (i * 1,2) are given in Appendix B. 

The homogeneous boundary conditions given by equations (21) are 
identically satisfied by the stress and displacement field equations 
(14-17). 

Following the same procedure given in the first section of this 

chapter the mixed boundary condition (22) can be written in the form of 

a singular integral equation given by 

h h c 

| (t^)" 2 dt + J Kn(x,t)v(t) dt + J Kj 2 (x,z) * 2 (z)dz 

-h -h 0 


up (x) ( 1 +k x ) , 

+ J K 13 (x,z)$ 3 (z) dz = - , | x I < h, (129) 

0 


where 


Kn(x,t) = J k n (x,t,n) e ^ t ^dn, (130. a) 

0 
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K i2 (x,z) = | k 12 (x,n) Cos x\z dn. 


(130. b) 


and 


Ki 3 (x,z) = | k 13 (x,n) Sinnz dn, (130. c) 

0 

kn(x,t,n) = { [ (a 22 -2a 21 ) Cosh nx - nxa 21 Sinh tjx] E 2 

+ [(-ai2 + 2ctn) Cosh nx + nxan Sinh tjx ] E 2 ]e ^ 

/(°n a 22"°i2 a 2i)) (131. a) 

ki 2 (x,il) = {“[(a 22 “2a 2 i) Cosh nx - nxa 21 Sinh nx] X 3 


+ [( - Oi2+2ctn) Cosh nx + nxan Sinh nx] X s }ne 

/(®ll®22”®12®2l) , 

k 13 (x,n) = {“[(a 22 -2o 21 ) Cosh nx - nxa 21 Sinh nx] \ k 

- [(-a 12 +2an) Cosh nx + nxo n Sinh nx] X s }ne 


-nh 


(131. b) 


-nh 


/(°i i<x 2 2 -Oi 2 a 2 1 ) , (131. c) 

where j , (i = 1,2, j = 1,2), and E^, (i = 1,2), are given in Appendix 
B. 

Similarly, by substituting equations (126) and (127) in equation 
(120) and separating the dominant part of the kernels, the following 
integral equation is obtained, 
h c 

J K 21 (y,t)v(t)dt + J K 22 (y ,z) $ 2 (z)dz - | $ 2 (y) 

-h 0 


+ J K 23 (y,z) $ 3 (z)dz = 0, 0 < y < c, 


(132) 


where 
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-i(y,t) = | kjj(t.n) e n ^ h t ^Cos ny dn, 


[(y.z) = | [k 22 (n) + 1] Cos ny Cos nz dn, 


K 23 (y,z) = | k 23 (n) Cos ny Sin nz dn, 


k 2l (t,n) = { [a 22 Cosh nh - nho 21 Sinh nh] 

+ [-o 12 Cosh nh + nhctji Sinh nh] E 2 }e ^ 

/(ouOjronOu), (134. a; 

k 22 (n) * {-[a 22 Cosh nh - nho 21 Sinh nh] X 3 

+ [-a l2 Cosh nh + nhaj ! Sinh nh] X 5 }ne ^ 

/ (<* i ia 2 2 ~Oi 2 a 2 1 ) , (134. b) 

k 23 (n) = {~l «22 Cosh nh - nha 2l Sinh nh] X* 

- [-a 12 Cosh nh + nhan Sinh nh] X s }ne"” I, k 

/(ai l a 22 -a 12 02 l ). (134. c] 

Similarly by substituting equations (126) and (127) in equation 
(123) and after separating the dominant part of the kernel, the follow- 
ing integral equation is obtained. 


| K 3 i(y,t)v(t)dt + | K 32 (y ,z) * 2 (z)dz + J K 33 (y,z)* 3 (z)d z 


- 2 ^a(y) = °> 0 < y < c, 
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where 


K a i(y,t) = J k 3 i(t,n) e" ,,(h " t) Sin ny dn, 
0 


(136. a) 


K 32 (y,z) = J k 32 (ri) 


(136. b) 


K 33 (y,z) = J [k. 3 3 (ti ) + l]Sin riz Sin tjy dn, (136. c) 

0 

k 3 X (y,t,n) = {[(o 22 -o 2 i) Sinh nh - a 21 T)h Cosh nh] E x 

+ [-Ou+aii) Sinh nh + a n i|h Cosh T)h] E 2 }e 

/(°i x a 22 ~a x 2 a 2 i ) , (137. a) 

k 32 (n) = {-[(o 22 -a 21 ) Sinh qh - a 21 nh Cosh i)h] X 3 

— Tl 

+ [(-ai 2 +a 11 ) Sinh nh + a lx nh Cosh T|h] X s } ne 

/(aua 22 -a 12 a 21 ), (137. b) 

k 33 (n) = {-[(a 22 -a 21 ) Sinh Tjh - a 2X nh Cosh i)h] X* 

- [(-ai 2 +a n ) Sinh uh + a lx nh Cosh t)h] X 6 }ne 


/(a 1 ia 22 -a X2 a 2 i), (137. c) 

and the functions a^j, (i = 1,2, j = 1,2), and E^, (i = 1,2), are given 
in Appendix B. 

Using integration by parts, the singular integral equation (129) 
can be obtained in terms of the slope of the crack opening displacement 
integral function 

GOO = (36) 


such that 
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h h c 

v(h) F(x) + J dt + j G(t) Kn(t,x)dt + J K 12 (x,z) * 3 (z) dz 

-h -h 0 

° f up(x)(l+Ki) 

+ J K l3 (x,z) d 3 (z)dz j— , I x I < h, (138) 

0 

Assume G(t) is a bounded function, then 

G(t) = H(t) , (139) 

where H(t) satisfies the Holder conditions in the closed interval 
[-h,h]. 

For c * 0, equation (138) reduces to equation (33) obtained for the 
case of the crack up to the interface. Hence the first three terms in 
equation (138) contribute the unbounded terms given by the left hand 
side of equation (80) expressed as 

Ui (x) = v(h) [B s e B ‘ (h-x) log(h-x) + B s e B ‘ (h+x) log(h+x) ] . (140) 

The unbounded terms contributed by the fifth term of the singular inte- 
gral equation (138) are given by 
c 

U*(x) = J kT 3 (x,z) ♦a(z) dz, (141) 

0 

where 

m 

kT 3 (x,z) = | k7 3 (x,n) Sin n z dn, (142) 

0 

and ki 3 (x,t|) is the asymptotic value of k 13 (x,i|) as i| + • and is given 

by 

k7 3 (x,n) = h(h-x) + Ni]e‘ n(h ‘ x) + [-n(h+x)+N 2 ]e _,,(h+x) 


(143) 


where N^, (i = 1,2), are given in Appendix B. Using the following 
results [23] 
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j I ,V'< h - x) Sin nz d, = 4L 


(144) 


K7a(x,z) = (h-x) 


2(h-x)z 

[(h-x) 2 +z 2 ] 2 


+ N 


l [(h-x) 2 +z 2 ] 


- (h+x) 


2(h+x)z 
[ (h+x) 2 +z 2 ] 2 


+ N2 [(h+x) 2 +z 2 ]‘ 


(145) 


Note that K i3 (x,z) becomes unbounded only as x ■+ ±h and z -*• 0 
simultaneously. 

Consider the following integral 
c 


I 


= J 


z $a( z) 
(h+x) 2 +z 5 


dz 


1 c . 0 3 (z) c * 3 (z) 

” 2 U z+i(h+x) + J z-i(h+x) ’ 

0 0 

According to [20, Chapter 4] near x = -h, 

I = -$ 3 (0) log(h+x) + bounded function. 

Similarly 

C . z * 3 (z) 

J dz = -* 3 (0) log(h-x) + bounded function. 

0 


(146) 


(147) 


(148) 


From equations (145), (147) and (148), the unbounded part contributed by 
the fifth term in the integral equation (122) is 

U 2 (x) = -0 3 (O) log(h-x) - $ 3 (0) log(h+x). (149) 
Following the same procedure, the fourth term in equation (138) is found 
to have no unbounded terms . 
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Adding equations (140) and (149), and using equations (123) and 
(124), 

Ux(x) + U 2 (x) = 0, (150) 

from which it follows that the sum total of the unbounded terms cancels 
out. This then proves that the assumption of a bounded slope function, 
G(t) is correct. 

Similar behavior of the slope function was reported by Erdogan and 
Backioglu [9] for the problem of a semi-infinite strip of finite width 
'2h' under a self -equilibrating pressure on the transverse end (y = 0) 
and stress-free ends (x = ±h). 

The axial stress in the half-plane has a logarithmic singularity at 
x = h, y = c. This is the well known logarithmic singularity [19, 24] 
due to the finite discontinuity in the shear load on a half-plane and is 
given by 

A 

lim a 2 (x,c) = K lim log (x-h) + Order (1), (151) 

x-h + 77 x-*h + 


A 

where K is the modified stress intensity factor and is defined by 


A - 9 y 3 

K = — o^y(h,c+) = — — [v 2 (h,c) - Vi(h,c)] 


4Vl 

= it (k x + 1) ♦ 3(c_) - 


(152) 


The expression (151) can also be derived by an asymptotic analysis simi- 
lar to that given in Appendix C. The axial stress in the half-plane at 
x = h, y = 0 is bounded, as the shear load is continuous and zero at 
that point. 
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The axial stress in the layer is also logarithmically singular at 

the interface crack-tip (x = h, y * c), and is given by 

A 

lim o l (x,c) = - K lim log (h-x) + Order (1). (153) 

x+h" 77 x-tr 

A 

Note that for a split length approaching zero, the modified SIF, K, is 
exactly half the value of the modified SIF, K given by equation (103) 
for the case of the broken layer (a = h). This is because the disconti- 
nuity in the shear load on the half-plane reduces by one-half for an 
infinitesimal split length. 

Normalizing the dimensions of equations (129), (132) and (135) with 
respect to the half-layer width 'h' and the half split length 'c ( , 


t 

s = h ■ 

(154. a) 

X 

r = h ’ 

(154. b) 

W = - , 

(154. c) 

c * 


(154. d) 

v(t) = e 1 (s)h. 

(154. e) 

*j(z) = Q*(w)c, 

(154. f) 

*a(z) * 0a(w)c, 

(154. g) 

p(x) = S(r) , 

(154. h) 


the following integral equations are obtained in a non-dimensional form. 
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J 9i(s) 

J (s-r) 2 
-1 


1 1 

ds + J h 2 Kn(hr,hs) 0i(s)ds + J c 2 K 12 (hr,cw)0 2 (w)dw 
-1 0 


irS(r)(l+K 1 ) 

+ J c 2 K 13 (hr,cw) 0 3 (w)dw — , | r | < 1, (155) 

0 

1 1 

J h 2 K 21 (cq,hs)0!(s)ds + j c 2 K 22 (cq,cw)0 2 (w) dw - y 0 2 (q) 

-1 0 


1 


+ 


/ 


c 2 K 23 (cq, cw) 0 3 (w)dw = 


0, 


0 < q < 1, 


(156) 


1 1 

J h 2 K 3 j (cq,hs)0j (s)ds + J c 2 K 2 3 (cq,cw)0 2 (w)dw 
-1 0 

1 

+ | c 2 K 33 (cq,cw)0 3 (w)dw - y 03(q) = °» 0 < q < 1. (157) 

0 


To find the numerical solution to the three integral equations 
(155-157), the following form of the unknown functions is assumed 
K 

0i (s) * l A,P.(s), (158. a) 

k=0 


L 

0 2 (w) = l B T (w), (158. b) 

£=0 


8s(w) 


M 

l 

m=l 


C m T m ( W )> 

m m 


(158. c) 


where P^Cs) and T^(w) are Legendre and Chebyshev polynomials of the 


53 


first kind, respectively. The normalized crack-opening displacement 
function 0 x (s) is approximated by a series of Legendre polynomials 
because of the availability of a direct relationship (109) for inte- 
grands with strong singularities. The other normalized displacement 
functions 0 2 (w) and 0 3 (w) are approximated by a series of Chebyshev poly- 
nomials. It has been shown [28] that a series of Chebyshev polynomials 
converges more rapidly than any other series of Gegenbauer polynomials, 
and converges much more rapidly than power series. 

Using the Hadamard concept, expressed mathematically by equation 
(103) for Legendre polynomials, the series approximations (158) for the 
unknown functions 0^(s), (i = 1,2,3), and the even and odd symmetry of 
0 2 (w)and 0 3 (w), respectively, equations (155-157) can be re-written as 

K 

I i (k,r) + 

k=0 

nS(r)(l+K 1 ) 

= AVU ’ 

K 

l A^Ck.q) + 
k=0 

I q I < i, 

K 

I A^Z 3 % (k,q) + 
k=0 


L M 

l B £ Z 12 (£,r) + l C m Z 13 (m,r) 

JL=0 m=0 

I r | < 1, (159) 

L M 

Y B £ Z 2 2(£ ) q) + Y. C m Z23(m,q) = 0, 

1 = 0 m =0 


( 160 ) 

L M 

X B £ Z 32 (£»q) + X C m Z 33 (m,q) — 0, 

£=0 m=0 


I q I < 1, (161) 

where Z^ j , (i = 1,2,3, j = 1,2,3) are given in Appendix B. 

Note that since the displacement function 0 3 (w) is an odd function 
of 'w' and is non-zero at w = 0, the first term of the series 
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approximation of the function 0 3 (w) is assumed to be of the form 
C 0 = C 0 sign (w), 

where 

sign (w) = + 1 , if w > 0 , 

= - 1 , if w < 0 . 

The collocation method is used to solve for the unknowns 
®1 an< * Sn' T * 16 c ^°^ ce collocation points is distributed 
symmetrically on the interval (-1,1) for equations (159-161). More 
points are concentrated near the end points to improve the rate of con- 
vergence. The roots of the Legendre polynomials satisfy such a 
selection of collocation points. Assuming r^, q^ and p as the kl^* 1 , 
il t ^ 1 and ml^ zero of the Legendre polynomials, Pk + i(s), Pl+i(w) and 
Pfl+l(w), respectively, equations (159-161) can be re-written as 

K L M 

X AfcZi i (k, r ki) X BjZi 2 (<l»rki) ^ C m Z 13 (m,r kl ) 
k =0 1=0 m =0 

irS(r )(1+ Kl ) 

= 4 ^ , (kl = 0,1,...,K-1,K), (162) 

K L M 

X ^^2 1 X 2(^>q£j) + X C^Z23 (m,q^^) = 0, 

k -0 1=0 m =0 

(11 = 0,1,..., L-l, L), (163) 

K L M 

X A^ZaiC^q^) + X B fc Z 32 (£,qmi) + X ^ m ^ 33 ^ m, %ii^ = 

k =0 1=0 m =0 


(ml = 0,1,..., M-l, M). 


(164) 
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Equations (162-164) represent a [(K+L+M+3) x (K+L+M+3)] system of 
linear equations which can be simultaneously solved for the unknown 
coefficients A^> B^, C^. The unknown function 0^, (i=l,2,3), can hence 
be calculated by using equations (158). Stresses and displacements can 
be evaluated at any point by using equation (126). 

The modified stress intensity factor given by equation (152) can be 
numerically expressed as 


4VUC M 


K = 


y c . 

L- m 


art) “ 


(165) 


The numerical scheme to evaluate the kernels Kjj, (i = 1,2,3, 
j = 1,2,3), and the finite integrals in equations (162-164) is explained 
in Appendix E. Convergence of results is also discussed. 



CHAPTER III 


RESULTS AND DISCUSSION 

The results presented here are for the case of plane strain with a 
constant pressure, p 2 , on the crack surface and no loads at infinity. 
This is the solution denoted by Ojj in equation (2). The complete solu- 
tion, ©total can be obtained simply by adding the uniform strain solu- 
tion, Oj, which is given by constant stresses as 


3 - 

II 

(166. a) 

(1-\>i 2 )E 2 

yy Pl (1 -v 2 2 )E 1 * 

(166. b) 


Five material combinations, given in Table I, are used in the 
results and are abbreviated as MC I, MC II, etc. for convenience. 


Table I. Material combinations used. 
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The Young’s moduli Ei and E* in the material combination MC I are 
the properties in the 90° and 0° directions of a unidirectional 
T300/5208 Graphite/Epoxy laminate. MC I is the most significant 
material combination because a transverse notch is more likely to 
develop first in a 90° ply. This behavior was observed by Masters [1] 
in his experimental study of impact loading of laminates. Most of the 
sample results shall hence be discussed for this material combination. 

The interleaf properties 2 were supplied by Masters as used in his experi- 
mental study [1]. 

Crack within the Laver 

In this problem, a symmetric transverse crack of length '2a' 

(a/h < l) 3 is situated in the layer and is opened by a uniform pressure, 
Px. The crack-tip axial stress has a square root type singularity and 
the coefficient of this singularity is defined as the mode I stress 
intensity factor (SIF), K x given by equation (61). When the SIF reaches 
a critical value Kx = K lc , the crack will propagate. This critical 
value is called the fracture toughness and is assumed to be constant for 
a particular material. It should be pointed out that the SIF depends on 
the crack length and the applied load. This criterion can be directly 
used only for brittle materials and needs modification for ductile 

1. T300/5208 Graphite/Epoxy composite laminate is made of 
T300-graphite fibers (manufactured by Union Carbide) in Rigidite 
5208-epoxy resin (registered trademark of Narmco Materials). 

2. Private communication of Professor J. G. Goree, Mechanical 
Engineering Department, Clemson University, Clemson, SC 29634 with Mr. 
John E. Masters, Chemical Research Division, American Cyanamid Company, 
Stamford, CT 06904. 

3. a/h = Crack length to layer width ratio. 
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materials, where yielding may exist. The simplicity of the criterion, 
however, makes it a single parameter to predict crack growth and frac- 
ture, and is useful even for cracks where yielding is in the form of 
small plastic zones. 

Before presenting the detailed results, the adequacy of modeling 
the interleaves as distributed shear and tension springs will be dis- 
cussed. Gecit and Erdogan [10] solved the problem of embedded cracks in 
periodic buffer strips separated by adhesive layers and modeled the lat- 
ter both as springs and as a continuum. They found small differences in 
the results obtained for the two models if the interleaf moduli and 
thickness were smaller than the strips. For example, for a material 
combination close to MC III and MC IV, and a/h =0.9 and (t/h < 0.2) 4 , a 
difference of less than 3% was found in the SIFs. Results from the 

Ki 

present study for the normalized stress intensity factor, — are 

Piv' a 

plotted as a function of interleaf thickness in Figure 5 for MC I. The 
values of the normalized SIF corresponding to the special cases of an 
interleaf thickness of zero as well as that of the thickness approaching 
infinity are found to be in excellent agreement with those special cases 
found in the literature. For the interleaf thickness approaching infin- 
ity, which is equivalent to an interleaf moduli of zero, or therefore 
corresponding to stress-free longitudinal sides of the layer, the normal- 
ized SIFs reach asymptotic values identical to those obtained by Sneddon 
and Srivastav [3] . They solved the problem of a transverse crack in a 
finite width strip with unloaded longitudinal edges. The normalized SIF 
for the other special case of an interleaf thickness of being exactly 


4. t/h = Interleaf thickness to half-layer width ratio. 
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equal to zero are found to be identical to those obtained by Hilton and 
Sih [5] and Bogy [6], They solved the problem of a layer, with an 
internal crack normal to the interface, bonded between two half-planes . 
These limit cases and the results of Gecit and Erdogan [10] for thin 
interleaves give considerable assurance as to the usefulness of the 
results in predicting the behavior of the interleafed composite lami- 
nates with an embedded crack. 

Additional significant results may be obtained from Figure 5, where 
it is seen that the normalized SIF increases with an increase in the 
interleaf thickness. This behavior is true for all material combina- 
tions. It does not, however, imply that the global strength of the lami- 
nate is necessarily decreased by inter leafing laminates. For example, 
if initial damage is in a 90° ply, the presence of the interleaf may 
result in a complete failure of the ply at a lower load. The subsequent 
behavior of the laminate when this ply is fully broken, and the influ- 
ence of the interleaf on delamination, additional ply failure, and ulti- 
mate laminate strength may be higher than with no interleaf. This 
behavior is considered in the next two sections. 

Further, from Figure 5, the SIF (numerator of the normalized SIF) 
is found to decrease as the crack length is increased for small inter- 
leaf thicknesses. Since the SIF is directly related to the load 
carrying capacity of a material , a decreasing SIF implies a higher load 
required to continue the crack. The crack growth is hence considered to 
be stable. For example, for a/h =0.9 and t/h < 0.08, the crack growth 
is stable. This behavior of a stable crack for small interleaf thick- 
nesses holds only for MC I, that is, for the case of the crack in a 
layer which is weaker than the half-plane. 
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a 2 (h,0) 

The maximum cleavage stress, in the half-plane occurs in 

line with the crack-tip at the interface and is plotted in Figure 6 as a 
function of interleaf thickness for MC I. For small crack lengths 
(a/h < 0.8) this stress is an increasing function in the range of small 
interleaf thicknesses (or large interleaf moduli) but is relieved on 
further increase in the thickness . This behavior is not found for large 
crack lengths. For example, for a/h = 0.9, the stress decreases mono- 
tonically with an increase in the interleaf thickness. Interestingly, 
for material combination MC V, where the crack is in a stiffer material, 
the increase in the cleavage stress in the range of small interleaf 
thicknesses is much more substantial for all crack lengths as shown in 
Figure 7. In such cases the introduction of a thin interleaf can, in 
fact, cause higher cleavage stresses and perhaps assist in continuing 
the crack across the interface. 

Figure 8 shows the transverse stress distribution along the inter- 
face for MC I. The stress is maximum and is tensile at y = 0, and 
decreases rapidly to become compressive away from the crack. Figure 9 
shows that the maximum transverse stress (at y = 0) decreases as the 
interleaf thickness is increased for MC I. This behavior is found to be 
typical and is exhibited for all material combinations. Table II gives 
the maximum interleaf transverse stress for different material combina- 
tions for a particular value of a/h = 0.9.' The percentage decrease due 
to the introduction of the interleaf is given in parentheses . 

Figure 10 shows the shear stress distribution along the interface. 
This stress is identically equal to zero at y = 0 and is very small near 
y = 0, but it changes sign and reaches a maximum value at about half a 
crack length away from the crack. The maximum interleaf shear stress, 














Table II. Maximum interleaf transverse stress for different material 
combinations (a/h = 0.9). 
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1 _ — 1 

Material 

Maximum Interleaf Transverse 
°xx(h’0)/P‘ 

Stress , 

Lomomat ion 

t/h = 0.0 

0.1 

0.2 

0.4 

MC I 

1.6889 

1.0358(39) 

0.7997(52) 

0.5776(66) 

MC III 

1.2942 

0.8753(32) 

0.7015(46) 

0.5239(60) 

MC IV 

1.2942 

0.2222(83) 

0.1350(90) 

0.0783(94) 

MC V 

0.3025 

0.1272(58) 

0.0893(71) 

0.0581(81) 


Table III. Maximum interleaf shear stress for different material 
combinations (a/h = 0.9). 



Material 

Combination 

Maximum Interleaf Shear Stress, 
1 °xv (h,y) 1 max/Pi 

t/h = 0 

0.1 

0.2 

0.4 

MC I 

0.4845 

0.3008(38) 

0.2404(50) 

0.1800(63) 

MC III 

0.3184 

0.2363(26) 

0.1996(37) 

0.1571(51) 

MC IV 

0.3184 

0.0741(77) 

0.0472(85) 

0.0285(91) 

MC V 

0.1168 

0.0374(68) 

0.0282(76) 

0.0196(83) 
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long the Interface for Constant Crack Length (MC I) 
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like the maximum interleaf transverse stress, is a decreasing function 
of interleaf thickness and is plotted in Figure 11 for MC I. Table III 
shows the same behavior for all material combinations. The percentage 
decrease due to the introduction of the interleaf is given in paren- 
theses . 


Crack Up to the Interface 

The only difference between this and the previous case is that the 
symmetric crack is now extended up to the interface and represents a 
broken layer. Due to the interleaf being modeled as a spring, the axial 
stress in the half-plane for this case has a logarithmic singularity at 
the crack-tip. A modified stress intensity factor, K is defined by the 
coefficient of this singularity and is given by equation (103). 

Figure 12 shows that this modified stress intensity factor K de- 
creases as a function of increasing interleaf thickness. This decrease 
is rapid up to t/h * 0.3. The modified stress intensity factor, K, un- 
like the stress intensity factor, Ki for the crack within the layer, 
does not directly give a measure of the load carrying capacity from the 
linear elastic fracture mechanics point of view. It may be recalled 
that modeling the adhesive as a spring resulted in a logarithmically 
singular axial stress in the half-plane at the crack-tip. A continuum 

model would give a power singularity but also not of a square-root 
type [10]. One of the methods to determine the strength of the compos- 
ite in such cases is to use failure criteria based on average stresses 
or point stresses and compare them with the ultimate strength of the 
undamaged material. A measure of the axial stresses at a specific point 
x = l.Olh in the half-plane as well as the average stresses calculated 


NORMALIZED INTERLEAF SHEAR STRESS 



NORMALIZED INTERLEAF THICKNESS, t/h 


Figure 11. Maximum Interleaf Shear Stress as a Function of Interleaf 
Thickness for Constant Crack Length (MC I) 



NORMALIZED MODIFIED STRESS INTENSITY FACTOR, K/p, 
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Figure 12. Modified Stress Intensity Factor as a Function of Interleaf 
Thickness for Different Material Combinations 
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over the interval x = h to l.Olh in the half-plane showed the same behav- 
ior of a rapid decrease as the interleaf thickness is increased to t/h = 
0.3. From the above observations and noting that the interleaves have 
an insignificant load carrying capacity and add to the weight of the 
structure, an interleaf thickness of approximately fifteen percent of 
the layer thickness seems suitable. It may be mentioned that graphite, 
epoxy and the interleaf have densities of the same order (0-5 pounds/ 
cubic inch) . 


H-Shaped Crack 

In this case, a symmetric transverse crack extending up to the 
interface and symmetric delamination along the interface is considered 
as shown in Figure 2. The axial stresses in the half-plane and the 
layer are, for this case, logarithmically singular at the interface 
crack-tips and are given by equations (151) and (153). Since the inter- 
leaf is modeled as springs, the transverse and shear stresses are finite 
in the interleaf even at the crack-tip. 

Figure 13 plots the transverse and shear stresses at the interface 
crack-tip, (x = h, y = c) as a function of split length for MC I. The 
transverse stress is tensile for no split and decreases rapidly as the 
splitting is initiated and grows. This peel (transverse) stress becomes 
compressive and therefore closes the split for relatively small split 
lengths. On the other hand, the interleaf shear stress at the interface 
crack-tip shows a very slow continued increase with split length growth. 
These results indicate that the split tip closes early and that most of 
the interface crack growth is due to shear alone. It is also found that 
the split length at which the crack closes decreases as the interleaf 
thickness is increased. The same behavior is exhibited for all material 
combinations. However, the order of the split length for the split to 



NORMALIZED INTERFACE CRACK TIP STRESSES 
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NORMALIZED SPLIT LENGTH , c/h 


Figure 13. Transverse and Shear Stresses at the Interface Crack-Tip as 
a Function of Split Length for Constant Interleaf Thickness 
(MC I) 
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close varies considerably for each material combination and is illus- 
trated in Table IV. The results are also compared with those obtained 
by the authors using the formulation and computer code developed by 
Goree and Venezia [12] for the case of two bonded half-planes with a 
transverse crack and delamination (T-shaped crack) along the interface. 

A further discussion on this behavior for longitudinal splitting is pre- 
sented by Wolla and Goree [29]. 

Table IV. Length of delamination until closure of crack for different 
material combinations. 


1 

Material 

Combination 

Split Length Until Closure of Crack, 
c/h 


— 

BBBB 

Results from f 12 1 

MC I 

0.35 

0.21 

0.260 

MC II 

0.018 

0.004 

0.096 

MC III 

0.205 

0.13 

0.164 

MC IV 


<0.002 

0.164 


The above results indicate that for a crack in the strong plies (MC 
II, MC IV), delamination (or at least positive peel stresses) is 
suppressed very early after initiation. For a crack in the weak plies 
(MC I, MC III), the extent of delamination is significant and suggests 
the need for the interleaves as well as that the bond-to-ply interfaces 
be of high tensile and shear strength. Also, from Figure 13 for zero 
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split length, the interleaf transverse and shear stresses at the inter- 
face crack-tip decrease with an increase in the interleaf thickness. 

This implies that a thicker interleaf will increase the external load 
required to initiate delamination. The same behavior is exhibited for 
all material combinations. 

Comparison with Experimental Studies 
Although results from a detailed experimental investigation are not 
available for direct comparison with the analytically predicted results, 
the preliminary experimental studies done by Masters [1] and Sun [2] are 
helpful in comparing the behavior of composite laminates with inter- 
leaves. These results support the analytically predicted behavior. Sun 
[2] used a baseline laminate specimen AS4/3501-6 Graphite/Epoxy 
(0 S /90 B /0 S ) and placed 5 mil adhesive film (FM1000 by American Cyanamid) 
between the 0° and 90° plies. The impact velocity required to initiate 
delamination in the adhesive layered (interleaf) laminate was twice as 
large as that required for the baseline laminate. Similar results were 
reported by Masters [1] for the AS4/1808 Graphite/Epoxy laminate 
[ (±45/0/90/0/90) 2 /±45/0/90/±45/]g. Figure 14 in reference [1], and 
reproduced here as Figure 14, shows the photomicrographs of the trans- 
verse cracks developed in the 90° plies of a baseline and an interleafed 
laminate. For the same impact loads, the interleafed laminate showed no 
delamination and damage was limited to transverse cracks. Delaminations 
occurred only for higher loads but were smaller than the ones developed 
for baseline laminates at comparable energy levels. 


ORIGINAL PAGE IS 
OF POOR QUALITY 


Figure 14. 
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BASELINE LAMINATE (200 X) 



INTERLAYERED LAMINATE (200 X) 

Transverse Cracks in Baseline and Interlayered AS4/1808 
Laminates at an Impact Energy Level of 800 in-lb/in 
(Reproduced from [1] by permission of Mr. John E. Masters, 
American Cyanamid Company, Stamford, CT 06904) 


CHAPTER IV 


CONCLUSIONS 


An analytical study is carried out to assist in the understanding 
of the influence of interleaves on the damage tolerance of multi-layered 
composite laminates. The geometry of the problem is idealized as a 
damaged layer bonded to two half -planes and separated by thin inter* 
leaves. The interleaves are modeled as distributed tension and shear 
springs. The damage in the layer is a symmetric transverse crack, which 
may extend up to the interface. Delamination along the interface is 
also analyzed. Material combinations assumed to approximate Graphite/ 
Epoxy laminates interlayered with thin thermoplastic film are used to 
discuss the results. The following conclusions are drawn from the study. 

A. For the case of a crack within the layer: 

1) The stresses at the crack-tip are singular and have, as 
expected, the classical square root singularity. 

2) The introduction of the low modulus interleaves increases 
the potential for the crack to extend while it is within 
the layer but reduces the stresses at the interface, 
which improves delamination damage tolerance. 

3) If the crack is in a layer which is stiffer than the 
half-planes, the use of interleaves results in higher 
cleavage stresses in the half-plane and may assist in 
continuing the crack across the interface. 

B. For the case of a crack up to the interface: 

4) By modeling the interleaves as distributed tension and 
shear springs, the axial stress in the half-planes is 
logarithmically singular at the crack-tip. 

5) The stresses in the interleaves and the half-planes 
reduce with an increase in the interleaf thickness. The 
rate of this reduction is most rapid for small interleaf 
thicknesses. For example, for the material combinations 
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of interleafed Graphite/Epoxy used in the results, this 
thickness is of the order of fifteen percent of the layer 
thickness. Any further increase in the thickness of the 
interleaf is an unprofitable addition to structural 
weight . 

C. For the case of an H-shaped crack: 

6) The axial stresses in the layer and the half-planes at 
the interface crack-tip are logarithmically singular, 
whereas the axial stresses in the half-planes at the 
intersecting points of the interface cracks and the 
transverse crack are bounded. 

7) The growth of delamination along the interface is stable. 
The split tip closes as the transverse stresses become 
compressive for small split lengths and further growth is 
due to shear alone. 

8) The split length until the split tip closes decreases 
with an increase in the interleaf thickness (or decrease 
in interleaf moduli). 

Based on the above conclusions, an interleaf thickness of the order 
of fifteen percent of the layer width is recommended for typical 
Graphite/Epoxy laminates. Since transverse cracks are more likely to 
occur in weak plies, the weak-stiff ply interfaces are prone to high 
stresses, which makes interleafing such interfaces a first choice. 
Similarly, the weak-weak ply interfaces are also prone to higher 
stresses but the possibility of a crossover of the transverse crack is 
also high. This makes interleafing such interfaces a secondary choice. 
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Appendix A 

Integral Expressions and Limit Evaluations 


The 
are used 


following integral formulas and limit evaluations 
in the formulation: 


f e Cos(ys)dy = 

dr“ J Q dr m 

oo 

,ro m 

— J e' Iy Sin(ys)dy = — 


[ r * + s 2 ^’ 


tr 2 + s 2 !’ 


r > 0, 


r > 0, 


oo 

J 

0 

00 

1 


1 

r(r 2 + s 2 ) 


r 

(r 2 + s 2 ) 2 


Sin(ry)dr = 777 (1 - e Sy ) , 


Sin(ry)dr = ~ y e 


I p “ Cos(ry)dr = ^ e‘ sy , 


I (r 2 + S 2 ) 2 Cos ( r y) dr = 4 ^ d + sy)e' sy , 
0 


J ( r 2 + s 2 ) 2 Cos ( r y) dr ■ O - s y)e" sy , 
0 


J e ry Cos rx Sin rt dr = *[( t y 2 + 

0 


t + x 
(t + x) 2 + 


, 27, 30] 


2] , y > 0, 
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j e' ry Cos rx Cos rt dr = i[ (t T + (t + x) 2 + y 2 ^ y > °» 


f e’ ry Sin rx Sin rt dr = ~ y 2 “ (t + x) 2 + ’ y > °* 

0 

Let a function f(t) be continuous on L and satisfy the Holder conditions 
on L + L ' , then 

oe 

iiS I f (t) !*“+ th ,)*] 1 * *t rt +T)’i dt = " ■ 


f fft'i f -* 3 - - y) 2 + -X 3 - x.(t ± Xl 2 .., 1 dt _ q 

J f(t) *[x 2 + (t - y) 2 ] 2 1 [x 2 + (t + y) 2 ] 2 * dt °» 
0 

oo 

f frt N f x 2 (t. + ,y) _ x 2 (t - y) , dt = 0 

J f(t) *[x 2 + (t + y) 2 ] 2 [x 2 + (t - y) 2 ] 2 * dt °* 


everywhere except at the end points of L. 
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Appendix B 

Definition of Constants 


The functions a^j, (i = 1,...,4, j 
(i = in equation (24) are 

an = - Cosh T)h, 
a J2 = ~ nh Sinh Tih, 

^ -Tih 

a 13 “ a l 4 — ® , 

Ml 

a 2 i = Sinh Tih, 

a 22 = Sinh Tih + uh Cosh qh , 


= 1» • • • >4) 


Mz _ 

a ” = uT * 


nh 


Mz -*h 

a 2 * = (Tih - 1) — e 1111 , 
Mi 


a 3 1 = - Sinh Tih, 

Kl - 1 


a 32 - 


2 

2m 2 t 


a 33 - - 


Sinh nh - nh Cosh nh, 

n - l, 


a 34 — 


2Mz t k 2 - 1 

[- -r- ^ - — r~ - e _T)h , 


a* i - 


-Cosh nh, 

Kl + 1 


a * 2 — 


Cosh nh - nh Sinh nh. 


a«,3 - [ 


a** — [■ 


2mz t n 

Ma 
2Mz t n 

Ms 


-nh 


k 2 + 1 


+ nh] e 


-nh 


, and D^, 


(nh - 1) - 


2 
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4 n* f ^ 

D,CJ>) = -J- J ( ^2 + pjl Cos £h d$, 
0 


4n 7 5* ♦»«) 

D 2 (ri) = “J J ( ^2 7 ^7j2 Sin $h d£, 
0 


2 7 * i ( 5 ) Ki - 3 

D 3 0 l) = -J J fp 7~^Tp [( — ) (£ 2 + Tl 2 ) + 2 U 2 ] Sin $h d$ , 

0 


2t| 2 . 0 l (£) Kl + 1 

D *<i) * * j Z(P + Py [( — — ) a 2 + tl 2 ) + 2 £ 2 ] Cos £h d 5 . 

0 

The expressions for the functions 04 j , (i = 1 , 2 , j = 1 , 2 ) in equa- 
tion ( 25 ) are 


on = [f CX 3 + X*) + 1] + [|(X, - X*)] e" 2l,h , 


X*i) 


“12 


n 2 h *'■ 1 

= [ 2 (Xs + x *) + ~T + ’ih + h 


2 h X*ri Xi 

+ [ V 1 ( " X ® + x *) ~ ~T + 2 “] 


-2nh 


2 2 


°2i = ^ ^*)] ^ [^("^s — X s ) + X 2 ] e ^ , 


n 2 i_ X *^ 1 

“22 = [ 2 (_Xs + X «) + ~ + 2] 


n 2u A 6 n A 2 « . 

+ [ 2 (X * + Xg) ~ x ** h “ ~r + 2~] e 2n > 


where the material constants Xi, (i = 1,...,6), are expressed as 
Klh2 “ K2hl 

X 1 = 


V2 + K 2 yj » 
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X 2 

X3 

X* 

Xs 

x« 


Pa “ Pi 
Pi + KiPa’ 


2PiPa 


K n (KaPi + 

Pa) 

2PiPa 


K s (R 2 yi + 

Pa) 

2PtP a 


K n (RiP 2 + 

Pi) 

2piPa 



K s (K,y 2 + Pi)' 


The above material constants are related by 
X 3 (X 2 + 1) = X«(X x + 1), 

XjXj = Xj,X s , 

which implies that there are four independent composite parameters. The 
dependence of the solution can hence be studied as a function of only 
four composite parameters. This concept has been discussed in detail by 
Dundurs [31]. 

The expressions for the functions F^, (i = 1,2), in equation (29) 

are 


Fi(n) = “X 3 t|Di + X*nD 2 - 


KaPi 

RaPi + Pa 


(Di ~ D 2 ) - 


Pa 

KaPi + y 2 


(D 3 + D„), 


™ = XsnDi + x 6 tid 2 + — ~ - (Di + d 2 ) + — - (d 3 - d„). 

The expressions for the functions H^n.t), (i = 1,...,4), in equa- 
tion (28) are 

Hi(n,t) * u[l - n(h - t) ] , 

H 2 (ri,t) = r\[Z - n(h - t)], 


84 


Ci - 3 

Ha (n » t) = n [T» (h - t) + — — ], 

K i +3 

H*(n,t) = n[-n(h - t) + — - — ]. 

The integrands E^, (i * 1,2), in equation (29) are 
EiCn.t) = l l [n 3 (h - t) (X 3 - X*) + T! 2 ( 2 X* - Xa) - X^], 

E 2 (n,t) = — f^i 3 0i - t) (-X s - X s ) + u 2 (X s + 2X S ) 

+ 2n 2 (h - t) Xj - 3X 2 n] . 

The functions f M and f 22 in equation (38) are 

fn = [-Xi + x*n] + I(-2Xa + 2X*)n 2 h - x*n + Xj] e" 2Tlh , 

f 22 = [-X 2 + x 6 n ] + [ ( 2 X S + 2X s )u*h - 4X 2 nh - x s n + X 2 ] e' 2n 

The material constants in equations (65) and (69) are 
Bi = C 2 + 4C g , 

B 2 = C s - 2C* , 

Ba = 2C 3 + 4C 7 - C 2 C« - 4C t 2 , 

B* = C% " 6C( , 

B« - C, - 4C S , 

B 6 = -C, + ✓ C, 2 - 4C 7 , 

where 

C 2 = T-— (“X B - X« + X 2 Xa + X 2 X*), 

A 3 A t 

2X 2 

3 X 3 X 6 * 

C* = r— — (-X 3 + X* + 3X S 4* 3X S - 4X 2 X 3 - 2\ 2 \*)> 

X 3 a 6 
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Cs — \ 3 x s (”XtX 5 + X 2 X t + 2 Xj - X 2 X 3 - X 2 X*), 

c. = (X. + x* + X s + X s ), 

c? = x7x7’ 

c, \ \ (X| * 2XjXg + X^Xg + 3X S • 2X2X3 X 2 X fc ) . 

a 3 a s 

The right hand sides of equations (126) are given by 
I*i (n) = (n) , 


I*a (n) — D 2 (n) , 

. • 

2y»n 

l* 3 (n) ■ D 3 (H) + — ^ J 0 2 (y) Cos ny dy, 

0 

2yiH ~ 

L*(n) = D*(n) + K ( + x J $ 3 (y) Sin ny dy, 

0 

and in terms of finite integrals as 


h 

Ui + 1) Li(n) = J v(t) Hi(n,t) e" l,(h-t) dt, 
-h 
h 

(Ki + 1) L 2 (n) = | v(t) H 2 (n,t) e _I * (h_t) dt, 

-h 
h 

( Kl + l) L 3 (n) = J v(t) H 3 (n,t) e -T|(h_t) dt 
-h 


2yin c 

+ — jj J * 2 (y) Cos ny dy, 

0 

h 

(Ki + 1) L*(n) = J v(t) H*(n,t) dt 

-h 


2yin c . 

+ x J *a(y) Sin ny dy, 

0 



The expressions for N^, (i = 1,2), in equation (143) are 


(X t ■ I 3 ■ X* * Xg) (h ~ x) 
N ‘ - ' 2 + * 

(X, - x 3 - x„ - x s ) (h + x) 


N 2 = -2 + 


2 X,X 


3*f 


The expressions for Z-y, (i = 1,2,3, j = 1,2,3), in equations 


(159-161) are 


Zn(k,r) = ^ {r Q k (r) - Q k+ i(r)} 

1 

+ | h 1 2 K lx (hs, hr) P k (s) ds, 

-1 


1 

Z X 2 (fc,r) = i J c 2 K x 2 (hr , cw) T fc (w) dw, 
-1 


(k=0, 1 , . . . ,K-1,K) , 


(fc=0 , 1 , . . . ,L-1 ,L) , 


1 

Z x 3 (0,r) = i | c J K x 3 (hr , cw) sign(w) dw, 
-1 


1 

Z X 3 (m,r) = i J c 2 K X 3 (hr, cw) T m (w) dw, (m=l, 2 , . . . ,M- 1 ,M) , 
-1 


1 

Z 2 x (k,q) = J h 2 K 2 X (cq, hs) P k (s) ds, (k= 0 ,l, . . . ,K- 1 ,K) , 

-1 

1 

Z 22 U,q) = | c 2 K 2 2 (cq, cw) T^w) dw - nc] , (fc= 0 , 1 , . . . ,L- 1 , 

-1 

1 

Z 23 ( 0 ,q) = i J c 2 K 23 (cq, cw) sign(w) dw, 

-1 


1 

Z 23 (m,q) = i J c 2 K 23 (cq, cw) T m (w) dw, (m=l, 2 , . . . ,M- 1 ,M) , 
-1 
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1 

Z 3 i(k,q) = | h 2 K 3 i(cq, hs) P k (s) ds, (k=0, 1, . . . ,K-1,K) , 

-1 

1 

Z 32 (£,q) = i J c 2 K 32 (cq, cw) T fc (w) dw, (£=0,1, . . . ,L-1,L), 
-1 
1 

Z 3 aC0,q) = «/ c 2 K 33 (cq, cw) sign(w) dw - uc sign(q)], 

-1 


1 

Z 33 (m,q) = i[ J c 2 K 33 (cq, cw) T m (w) dw - ire], (m=l,2, . . . ,M-1,M) . 
-1 


Appendix C 

Stress Intensity Factor for the Embedded Crack 
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The stress intensity factor for an embedded crack (a < h) is given 
by equation (61) as 


lim 


K x = / 2(x - a) a 1 (x,0) . 

1 x-ka^ yy 5 ' 

From equations (30-41), the axial stress q ^ (x,0) can be written as 


(61) 


o 1 (x,0) 

yy 


fGOl 

1 + Ki J t - > 


dt + bounded function. 


-a 


G(t) has an integrable singularity at t = ±a given by 
H(t) 


G(t) = 


/ a 2 - t 2 


or 


G(t) = 


Hft) e 


■a i/2 


(t 


- a)^ (t + a)^ 


Consider the sectionally holomorphic function 

Gill 


«»> 


-a 


(C.l) 


(51) 


(C.2) 


(C.3) 


According to Muskhelishvili [20, Chapter 4] and using equation (C*2) 

i n/2 


0(z) = 


H(-a~) e 


H £a) 


/2a / z + a /2a /z-a 


+ $o(z), 


(C.4) 


where $ 0 (z) is a bounded everywhere except at the end points z = ±a, 
where 

H 0 (±a) 


z ± a 


*o 


, Re do) < i. 


*o(z) | £ 


(C.5) 
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Using the Plemelj formulas (47) on equation (C,4), the function 
0(x) may be defined near x = a by 


0(x) = - 


H(a) 


+ *o(x). 


/2a / x - a 

where $ 0 (x) is a bounded function and has a behavior similar to 0„(z) 
at the end points. From equations (C.l), (C.3), (C.6), and (61), the 
value of the mode I stress intensity factor is given by 
4 Hi 


(C.6) 


Ki = - 


H£jQ 

(1 + Ki) /a 


a < h. 


(C.7) 


In terms of the non-dimensional function ’f(s), defined by equation (56), 
the stress intensity factor, K x , can be re-written as 


Ki 


4Vl Ul) 
(1 + ) / a 


a < h. 


(C.8) 



90 


Appendix D 

Cleavage Stress in the Half-Plane 
for the Broken Laver 


The cleavage stress in the half-plane is given by 
4p 2 “ r 

o£y(x,0) = — J [fi(n) + (UX - 2)g,(n)] e _Tlx dn. (D.i) 

0 

From equation (24), the functions f 2 (n) and g 2 (n) may be expressed as 
— f 2 (n) e 1,11 = f i (n ) [-T)h e *** 1 + Cosh nh] + gi(t|) [~n*h 2 e* 1 * 1 ] 

Ul 

+ (1 - Hh) Dx + nh D 2 , CD. 2) 


“ g 2 (n) e~ T ' h = fi(n) [e 1,h ] + gx(n) [ "nh e^ + Sinh nh] 

+ D t - D 2 . (D. 3) 

Substituting equations (D.2), (D.3), and (24) in equation (D.I) for 
a = h, 

4v ** h 

o^y(x, 0 ) = | j (SxTi 2 + 6 2 n + 6 3 ) v(t) e " n(x ~ t) dt dn 

0 -h 

4hx ~ 

+ — J {[-tl(h-x) - 1 ] Dx 
0 

+ [n(h-x) + 2] D 2 } e _Tl(x_h) dn 

+ bounded function, (D.4) 

where 

«x = Qi(h - t)(h - x) + (h - t) - (h - x), 

62 = -Q, Qx(h - t)(h - x) + (Q 3 - Q g )(h - t) 

+ (Q* -Qs)(h - x) - 3, 


91 


and 


63 = (-Q10 Qi + Qs 2 Qi)(h - t)(h - x) 
+ (Qs - Qio - Qs Qa + Q 9 2 )(h - t) 
+ (Qs + Qio - Qs Q* + Qs 2 ) (h - X) 
+ (Qt + 3 Qs ) • 

Ql “ » \ (** X 3 1" 1* " ^3^2 ^ 2 ^s)» 

A 3 A 6 


Q 3 “ o\ \ (~ 3 X 3 "*■ 3X4 + 2X2X4 — 4 X 2 X 3 + Xg + X 5 ) , 

ZA 3 Ag 


Q* 

Qs 

Qs 


2 x 3 x 6 + ^1^6 

_ A. 

X 3 X g * 

A 

x 3 x s ’ 


• 4X4 * 4 " 2X3 + 3X3X3 — 3X3X11), 


q 7 = 


2X3X5 


(2X5X5 + X5X5 “ 6X4 + 3X3 ■ 3X3X4 + 6X3X3 * X5 2X5), 


Qs = (Xi + x 2 ). 


Qs = 2)^x7 + X* + X s + X $ ), 


Qio — 


X 3 X 8 


The values of X^, (i = 1,...,6), in terms of material constants are 
given in Appendix B. 

Using the following relationships 




1 dx 


< h - *> k <;rh> - (h - *> h <;rh> - ;rh- 


C ' cX- 
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(h - « 5F «rh> ■ < h - « 5? <rh> - 4 c^. 


h 

r v(t') 

J t - x = v (h) (x - h) + bounded terms, for x ■+ h, x > h, 
-h 


find equation (D.4), the cleavage stress is given by 



4y 3 v(h) 

log (x - h) + Order (1). 


(D.5) 


93 


Appendix E 

Numerical Techniques and Convergence of Results 

All the computational work was done in double precision on the NAS 
AS/XL-60 computer at Clemson University. The computer has a double pre- 
cision unit of EPS — 0.222 x 10" 15 and can represent a maximum of 18 
decimal digits. All programs were written in FORTRAN. 

Crack within the Layer 

The Fredholm Kernel in equation (57) is an integral over a 
semi-infinite interval and is estimated by using the NAG [32] FORTRAN 
library routine called D01BAF. It computes an estimate of the integral 
using the n-point Gauss-Laguerre formula [25], which is exact for any 
function of the form 
. 2n-l 

f(x) = e“ X £ CjX 1 (E.l) 

i=0 

This routine is appropriate for exponentially decaying functions and the 
parameter 'b' is chosen to match the decay rate of the function, f(x). 

The stress intensity factor, K t given by equation (61) is the most 
significant parameter and is tabulated to 6 significant digits as a 
function of the number *N' of the Gauss -Chebyshev quadrature formula 
used in equation (57). Results are given for material combinations MC I 
and MC V in Table E-I (~ = 0.9, ~ = 0. 1) . It may be mentioned that the 
techniques [22] used to evaluate the stress intensity factor in this 
case are based on treating the singular integral equation as a system of 
Fredholm equations and then finding the unknown function at discrete 


collocation points. 
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Table E-I. Stress intensity factor as a function of number of quadra- 
ture points (N) . 


N 

Normalized Stress 
Kx/p 

Intensity Factor, 
x 'f a 

Material Combination 

MC I 

MC V 

25 

2.05261 

0.847155 

29 

2.05262 

0.847154 

33 

2.05262 

0.847154 

37 

2.05262 

0.847154 

41 

2.05262 

0.847154 


Convergence of the above results needs no explanation. Fortv-one 
quadrature points were used for all combinations. 

Crack Up to the Interface 

The kernel K u in equation (106) is solved by the techniques 
explained in the previous section. The finite integral in equation 
(106) given by 
1 

I = J h 2 P n (s) K ll (hs,hr i ) ds (E.2) 

-1 

is solved by using the n-point Gauss-Legendre formula [25] exact for any 
function of the form 
2n+l 

f(x) = l Ci x 1 . (E. 3) 

i=0 

The formula is appropriate for functions which can be approximated by 
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polynomials, and is given by 


b ' 

I f(x) dx 
a 



m 


l 

i=l 


Wi 


f(yi)» 


(E.4) 


where 


Yi « ( 


b -j, 


Xi + 


,b + a 


), 


Xi = ith zero of P m (x) , 

»i • (i . 2 X1) « «•= 

and P m (x) is a Legendre polynomial. 

The stress intensity factor, K is the slowest converging parameter 
in the problem and is tabulated in Table E-II as a function of the num- 
ber of collocation points, N for material combinations MC I and 
MC II (J = 0.2). 


Table E-II. Modified stress intensity factor as a function of number of 
collocation points (N) . 


N 

Normalized Modified Stress Intensity Factor, 

K/pi 

Material Combination 

MC I 

MC II 

24 

0.69068 

0.79443 

28 

0.68820 

0.79321 

32 

0.68650 

0.79238 

36 

0.68528 

0.79178 

40 

0.68437 

0.79132 
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Analytically more accurate results can be obtained by increasing 
the number of collocation points, N. This is limited only by computer 
time, but for practical purposes and within reasonable accuracy, a 
limiting value of N is required. 

One can estimate the value of the normalized stress intensity 
factor, K/px as N-+- by using extrapolation techniques [22,33] for 
accelerating the convergence of a monotonic sequence (oscillating 
sequences require modification) . The normalized stress intensity 
factor, K/px is expressed as a function of the number of quadrature 
points, N as 

K(N) = K(-) + S ~, (E.6) 

N° 

where K(~) = estimated value of K/px, 

B = unknown constant, 
a = exponential constant. 

The three unknowns can be evaluated by using three values of N (N^, 
i = 1,2,3). To ensure stability of the calculated results, the exponen- 
exponential constant, a should be greater than one. Some tabulated 
values up to four significant digits are given in Table E-III for 
material combination MC I (jj * 0.2). 

Table E-III. Estimated values of the normalized modified stress 
intensity factor. 
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For material combinations MC I - MC IV, the difference between the 
estimated value, K(-) and the calculated value, K(40) for forty colloca- 
tions points is less than one percent. Hence forty collocation points 
are used for all the solutions. For the material combination V, conver- 
gence is slower and requires more than forty collocation points to 
achieve the same accuracy. This difficulty can only be overcome at the 
expense of computer time. 


H-Shaped Crack 

The kernel K u and the finite integral in equation (162) are solved 
by using the techniques explained in the previous sections . The kernel 
K 12 given by equation (162) is evaluated by the following procedure for 
speed and accuracy. The kernel K 12 can be written as 

K 12 (x,z) = kT 2 (x,z) + K^ 2 (x,z) + K^ 2 (x,z) , (E.7) 

where 


k7 2 (x,z) 


Ki’ 2 (x,z) 


K^ 2 (x,z) 


J k7 2 (x,n) Cos nz dn, 
0 

J ki‘ 2 (x,ri) Cos nz dn, 
0 


f «o L 

j [ki 2 (x,n) - k j 2 (x,n) - ki 2 (x,n)] Cos nz dn. 

0 


(E.8a) 


(E.8b) 


(E. 8c) 


and kt 2 is the asymptotic value of the integrand ki 2 (x,n) as n^* and 
x ■* ±h and is given by 

k72(x,n) = [(din + d 2 ) Cosh nx + (d 3 n + d k ) Sinh nx] e T,h , (E.9) 
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and k^ 2 (x,ri) includes the terms of order (1/n) as n ■+ ~ and x -*• ±h, and is 
given by 


k ^ x > 1 » = [^ m ; + vttJ Cosh 1,31 e " nh 


d 9 


d 1( 


+ [ — — - + — — H Sinh i)x e T)h , 
L H + d e ri + d, J 


(E. 10) 


where d^, (i = are explicit expressions in terms of material 

properties and geometrical parameters . 

KT 2 can be evaluated exactly by using the following exact result 

[27], 


[ e Cos T)z dr) = T“i [ 2 , — 2 ] , 

, m J dx m l x + z* 1 ’ 

dx „ 


x > 0. 


(E.ll) 


K ^2 is an integral of a slowly converging function [0(~)] over a 
semi-infinite interval. The following results express the integral K^ 2 
as finite interval integrals and a highly convergent semi-infinite inte- 
gral . 

00 b 

„ -ax . -ax e . 

J ~ Cos cx dx = e ab Sin (be) [Tan - * ~ - J x in CX dx] 

0 0 

+ e ab Cos (be) [~i log ( C — ~T '~ ) 

a 


/ 


1 - Cos cx -ax 


e dx + Ei (ab) ] , 


a > 0, b > 0 

where E 2 (x) is the exponential integral defined by 


(E. 12) 


-u 

E t (x) = / du, x > 0. 


(E«13) 


x 
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The finite integrals in equation (E.12) are estimated by using a 
NAG [32] library routine D01AKF, an adaptive integrator suitable for 
oscillating integrands. It uses Gauss 30-point and Kronrod 61-point 
rules. Another NAG [32] library routine S13AAF calculates an approxi- 
mate value of the exponential integral Ei(x). 

The break-up of the integral as given by equation (E,7) allows the 
kernel K l2 to be evaluated accurately and rapidly because the integrand 
of K^ 2 decays at least at the rate of e T1 ^ 1 , where 'h' is the half layer 
width. K^ 2 is evaluated by using a NAG [32] library routine D01ANF, an 
adaptive integrator for finite cosine or sine transforms. In view of 
the fast convergence [0(e T1 ^ 1 )] of the integral and definition of K^ 2 > 
the upper limit in the integral can be substituted as 200/h. 

The rest of the kernels in equations (162-164) are evaluated 
similarly by breaking up the integrand in form of equation (E»7) and 
using the following additional results. 


,m . ,m 

— - I e " X Sin qz dq = — “ ("i T - ;*), x>0, (E.14) 

dx m J Q dx m X + 2 

«» b 

J ^ + k Sin cx dx = e ab Cos (be) [Tan”* ^ - J “ CX dx] 

0 0 

- e ab Sin(bc) log (- — ^T 9- ) 

b 

+ J e dx + Ei(ab)] , 

0 

a > 0, b > 0, (E.15) 

ee 

I IPTb ^ = ‘ Sin(ab > [ Si(ab > " f3 

o 


- Cos(ab) Ci(ab), a > 0, b > 0, 


(E* 16) 


I TV*b ^ = Sin(ab) Ci < ab ) 

o 

- Cos(ab) [Si(ab) - |], a > 0, b > 0. (E.17) 

The NAG [32] library routines S13ACF and S13ADF estimate the value 
of the cosine integral Ci(x) and the sine integral Si(x), respectively 
which are defined [25] as 

x 

Ci(x) = y + log x + | - gs * — - dt, x > 0, (E. 18) 

0 

X 

Si(x) = | t dt, x > 0, (E. 19) 

0 

where Z = Euler’s constant (= 0.577215 ...). 

Equations (E.12) and (E.15) may give erroneous results for | ab | > 5 

due to round-off errors. For such cases, the integral is directly 

solved by using the NAG [32] library routine, D01ANF, an adaptive inte- 

integrator for cosine or sine transforms of a function over a finite 

grator for cosine or sine transforms of a function over a finite inter- 
val. The upper limit in the integral is taken as 200/b, where b is 
the exponential decay rate of the integrand K 12 . 

All the finite integrals in equations (162-164) are solved by using 

the technique explained in the second section of this appendix. 

To illustrate the convergence of the results, the normalized stress 

A 

intensity factor, K/pi is tabulated as a function of the collocation 
points of the quadrature scheme used in equations (162-164) for 
material combinations MC I (^ = 0.05, - = 0.2) and MC II (“ = 0.002, 
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Table E-IV. Modified stress intensity factor as a function of number of 
collocation points (K, L, M) . 


K=L=M 

Normalized Modified Stress Intensity Factor, 

ft/px 

Material Combination 

MC I 

MC II 

24 

0.36473 

0.39367 

28 

0.36461 

0.39320 

32 

0.36454 

0.39296 

36 

0.36450 

0.39287 

40 

0.36447 

0.39287 


Applying the same extrapolation technique and accuracy criterion 
illustrated in the previous section, forty collocation points (K=L=M=40) 
are found to be sufficient for material combinations MC I - MC IV. For 
the material combination MC V, the convergence is slower and requires 
more than forty collocation points. This difficulty can be overcome at 
the expense of computer time. 










Appendix F 
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CYJkNJUHID 


American Cyanamid Company 
Stamford Research Laboratories 
1937 West Main Street 
P.O. Box 60 

Stamford, Conn. 06904 - 0060 
Tel. (203) 348-7331 


June 18, 1987 


Mr. Autar K. Kaw 

Department of Mechanical Engineering 
318 Riggs Hall 
Clemson University 

Clemson, South Carolina 29634-0921 
Dear Mr . Kaw : 

Please excuse my tardiness in issuing this formal, 
written response to your request for photographs. As I indicated 
during our telephone conversation, including these photographs 
(i.e. Figs. 2, 13 and 14 from my paper entitled "Characterization 
of Impact Damage in Graphite/ Epoxy Laminates") in your report 
presents no problem. I hope this letter will meet your 
requirements for a formal consenting letter. If not, please call 
and I will be happy to respond. 

As we also discussed during our^ telephone conversation, 
I would be most interested in receiving a copy of your final 
report . 


Sincerely yours, 



J EM ) saj 
Doc. (1.20) 
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